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We have studied the quantum oscillations of the conductance for arrays of connected mesoscopic 
metallic rings, in the presence of an external magnetic field. Several geometries have been considered: 
a linear array of rings connected with short or long wires compared to the phase coherence length, 
square networks and hollow cylinders. Compared to the well-known case of the isolated ring, we 
show that for connected rings, the winding of the Brownian trajectories around the rings is modified, 
leading to a different harmonics content of the quantum oscillations. We relate this harmonics 
content to the distribution of winding numbers. We consider the limits where coherence length 
Lip is small or large compared to the perimeter L of each ring constituting the network. In the 
latter case, the coherent diffusive trajectories explore a region larger than L, whence a network 
dependent harmonics content. Our analysis is based on the calculation of the spectral determinant 
of the diffusion equation for which we have a simple expression on any network. It is also based 
on the hypothesis that the time dependence of the dephasing between diffusive trajectories can be 
described by an exponential decay with a single characteristic time Tip (model A) . 

At low temperature, decoherence is limited by electron-electron interaction, and can be modelled 
in a one-electron picture by the fluctuating electric field created by other electrons (model B). It 
is described by a functional of the trajectories and thus the dependence on geometry is crucial. 
Expressions for the magnetoconductance oscillations are derived within this model and compared to 
the results of model A. It is shown that they involve several temperature-dependent length scales. 

PACS numbers: 73.23.-b ; 73.20.Fz ; 72.15.Rn 



I. INTRODUCTION 



Understanding which processes Hmit phase coherence 
in electronic transport is an important issue in meso- 
scopic physics. Such phenomena hke weak locahzation 
or universal conductance fluctuations are well understood 
to result from phase coherence effects limited at a given 
time (or length) scale called phase coherence time T,p (or 
phase coherent length L^p). This explains the interest in 
studying quantum corrections to the classical conductiv- 
ity : to provide a powerful experimental probe of phase 
coherence in weakly disordered metals and furnish some 
informations on the microscopic mechanisms responsible 
for the limitation of quantum coherence. This limitation 
originates from the coupling of electrons to external de- 
grees of freedom like magnetic impurities or phonons^'^. 
It also results from the interaction among electrons them- 
selves. The physical origin of this decoherence in weakly 
disordered metals has been understood in the pioneer- 
ing paper of Altshuler, Aronov & Khmelnitskii (AAK)P1. 
In a one electron picture, it is due to the fluctuations 
of the electric field created by the other electrons. In a 
quasi-ld wire, these authors have shown that this mech- 
anism leads to the following temperature dependence of 
the dephasing time t,^{T) oc T^^/^ . This power-law can 
be understood qualitatively as follows: the typical de- 
phasing is proportional to the fluctuations of the electric 
potential, which themselves are known from Nyquist the- 
orem to be proportional to the temperature T and to the 



resistance of the sample. For an infinite wire, the rele- 
vant fluctuations are limited to the scale of the coherence 
length itself. Consequently, the dephasing time has the 
structure : h/T^p = k^T / g{L^), where g{Lip) is the di- 
mensionless conductance at the length scale L.^. For a 
quasi-ld conductor, the conductance is linear in length 
and the length scales as the square-root of time. There- 
fore the function gi^L^p) scales as y/r^, whence the above 
power law. 

More recently, Ludwig & MirliiP and two of the 
authors^ have considered the geometry of a ring, and 
they have shown that the damping of magnetoresistance 
oscillations could be described with a different temper- 
ature dependence of the dephasing time T^piT) cx T~^. 
This new behaviour can be qualitatively understood by 
considering that the diffusive trajectories encircle the ring 
and have all a length equal to the perimeter L of the ring, 
so that the relevant resistance is the resistance of the ring 
itself. As a result : h/rp = kBT/g{L). 

In Ref.El we have shown how the dephasing on a ring 
depends on the nature of the diffusive trajectories : the 
fluctuations of the electric potential affect differently tra- 
jectories which encircle the ring and trajectories which do 
not encircle it. Within this framework, we have analyzed 
magnetoresistance experiments performed on a square 
network of quasi-ld wires, and we have found that in- 
deed two characteristic lengths with two different tem- 
perature dependence could be extracted from the datcP. 
These recent considerations have led us to the general 
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conclusion that the dephasing depends on the geometry 
of the system considered. 

The purpose of this paper is to analyze the dephasing 
process and to calculate the weak localization correction 
in different geometries, where the decoherence induced 
by electron-electron interaction may have a more com- 
plex structure. In order to address this question, it is 
important to understand that the weak localization cor- 
rection depends on two ingredients, one is the probability 
to have pairs of reversed trajectories, which is related to 
the return probability V(t) for a diffusive particle after a 
time t, the other is the nature of the dephasing process 
itself. Schematically, the weak localization correction to 
the conductivity can be written as 



dtV{t) (e 



i>(t)\ 



(1) 



where (e'*^*-') is the average dephasing accumulated 
along a diffusive trajectory for a time t. The return 
probability has been analyzed in Ref .'^ for various types 
of networks. Its Laplace transform, the spectral deter- 
minant, can be simply calculated from the parameters of 
the network. More complex is the analysis of the dephas- 
ing process itself. A simple and natural ansatz would 
be to assume an exponential decay (e'**^*-') = e~*^'^f. 
This assumption is correct when the dephasing is due 
to random magnetic impurities or electron-phonon scat- 
tering. For electron-electron interaction the analysis of 
the AAK result for a wire shows that time dependence 
is not exponentiaP : (e'**^*))ee G"'/'^'^ . The qualitative 
reason stands again on the fact that dephasing can be de- 
scribed as due to the fluctuations of the electric potential 
due to other electrons. Then, one may understand that 
this dephasing depends on the nature of the trajectories 
and is not exponential. The main goal of this paper is 
to describe this dephasing for complex networks and to 
generalize the known results of the infinite wire and the 
ring. 

The paper is organized as follows : In section |ll] we 
recall the physical basis at the origin of this work and in 
section |III| we present the general formalism appropriate 
for our study. In the next sections, we consider succes- 
sively more and more complex geometries. In section |IV] 
we recall known results for the infinite wire and the ring. 
In section |V] we consider the case of a ring attached 
to arms and show how the harmonics of the magnetore- 
sistance oscillations are reduced by the existence of the 
arms. The situation is the same for a chain of rings con- 
nected through arms longer than the coherence length. 
When rings become close to each other the dephasing in 
one ring is strongly modified by the winding trajectories 
in the neighboring rings. This is discussed in section [VT] 
The case of an infinite regular network is much more 
difficult to address since the hierarchy of diffusive trajec- 
tories is difficult to analyze, and we have used the limit 
of the infinite plane as a guideline (section VII I . Finally 
the case of a hollow cylinder (section |VIII ) is quite in- 



teresting since it combines trajectories winding around 
the axis of the cylinder and two-dimensional trajectories. 
Throughout the paper, we shall consider two situations, 
respectively denoted by model A and model B : the case 
where the dephasing has a simple exponential time de- 
pendence, and the case where the dephasing is induced 
by electron-electron interaction. We shall systematically 
discuss the analogies and the differences between these 
two situations. 



II. BACKGROUND 

In a weakly disordered metal, due to elastic scattering 
by impurities, the classical conductivity reaches a finite 
value at low temperature, given by the Drude conduc- 
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tivity ctq = "'^m" ' "^ti^re n is the electronic density and 
Te the elastic scattering time. Quantum interferences are 
responsible for small quantum corrections to the Drude 
result. One important contribution, that survives aver- 
aging over the disordeJ ^ * ^" * ^^ !, comes from interferences 
of reversed closed electronic trajectories, and therefore 
diminishes the conductivity. This quantum contribution 
to the average conductivity is called the weak localiza- 
tion (WL) correction. It has been expressed as (jlj where 
the function (e'* '■*•') describes dephasing and cut off the 
large time contributions. A simple exponential decay 
^gi$(t)^ _ Q~t/T^ jg usually assumed (denoted model A 
in the present paper), is the phase coherence time, re- 
lated to the phase coherence length L^p — yj Dr^p, where 
D is the diffusion constant of electrons in the disordered 
metal. From eq. (fT|, we obtain the WL correction in a 
wire Acr'^'^) ^ — J^*" ^ ^ ^ —L^ and in a plane 

Acr(2d) j^-f dl ln(T^/re) (diffusion sets in after 

a time Te, whence the lower cutoff in the integrals). In 
practice, the WL is a small correction to the Drude con- 
ductivity and it can be extracted thanks to its sensitivity 
to a magnetic field. In the presence of a magnetic field, 
the contribution of a closed diffusive trajectory C is mul- 
tiplied by e^"^'^'^/^, where 0c is the magnetic flux through 
the loop. This phase factor comes from the interference 
of the two reversed electronic trajectories, whence the 
factor 2. After summation over all loops, the additional 
magnetic phase is responsible for the vanishing of the con- 



tributions of loops such that 



> 



where 



h/e 



is the flux quantum. Therefore the magnetic field pro- 
vides an additional cutoff at time corresponding to 
diffusive trajectories encircling one flux quantum. In a 
narrow wire of width w we have r. 



(Id) 

B 



bl/{Dw^B^) 

and in a thin film (plane) t^'^^ ^ (f)o/{DB) (see RefeP^l), 
The two cutoffs are added "a la Matthiessen"'^2llSl as 
1 /r^ 1 /t^ + ^/ti3 ; this leads to a smooth dependence 
of the WL correction as a function of the magnetic field. 

The above discussion concerns homogeneous devices 
(like a wire or a plane). Another experimental setup ap- 
propriate to study quantum interferences and extract the 
phase coherence length is a metallic ring or an array of 
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rings. In this case the topology constrains the magnetic 
flux intercepted by the rings to be an integer multiple of 
the flux per ring (j) (we neglect for the moment the pene- 
tration of the magnetic field in the wires) : (f>c = ncj) with 
n e Z. This gives rise to Aharonov-Bohm (AB) oscilla- 
tions of the conductance as a function of the flux with 
period ^q. Disorder averaging is responsible for the van- 
ishing of these 0o-periodic oscillations : only survive the 
contributions of the reversed electronic trajectories lead- 
ing to WL correction oscillation s, kno wn as Al'tshuler- 
Aronov-Spivak (AAS) oscillation^iSIIll^ with a period half 
of the flux quantum. It will be convenient to introduce 
the harmonics Acr„ of the periodic WL correction. An 
important motivation for considering the harmonic con- 
tent Acr„ in networks, is that it allows to decouple the 
two effects of the magnetic fleld^^ : the rapid AAS os- 
cillations (A(T„^o) Sind the penetration of the magnetic 
fleld in the wires, responsible for a smooth decrease of 
the MC at large field (Actq). Since the n-th harmonic is 
given by contributions of loops encircling n fluxes we can 
write 



2e^D 



TTS 



dtVn{t) e 



-t/r 



(2) 



where Vn{t) is the return probability after a time t having 
encircled n fluxes, s is the cross section of the wires. In 
an isolated ring of perimeter L, this probability reads 



\/4TrDt 



exp - 



Act, 



AAS 



(nLT 
4Dt 



2e2 



Integral 



21 givePI 



Lip G 



-\n\L/L^ 



(3) 



where ~ \/Dt^ is the phase coherence length. Note 
that Act_„ = Act„ follows from the symmetry under re- 
versing the magnetic field ; in the following we will simply 
consider n ^ 0. The exponential decay of the harmon- 
ics directly originates from the diffusive nature of the 
winding around the ring : for a time t t^, the typi- 
cal winding scales as nt ~ y/Dt/L ~ L^/L. The AAS 
oscillati ons we re first observed in narrow m etallic h ollow 
cylinderPE^ and in large metallic networkPE^Iini 

Although the simple behaviour (|3| has been used to 
analyze AAS or AB oscillations^^ i n ma ny experiments 
until recently (see for example Refs.li^Ml)^ a realistic de- 
scription of a network made of connected rings leads to 
harmonics with a LjL^ dependence a 'priori quite dif- 
ferent from the simple exponential prediction (|3| for two 
reasons related to the nontrivial topology of the networks. 

(i) Winding properties of diffusive loops in networks. - 
Consider for example the square network of figure [l] made 
of rings of perimeter L = 4a. For <C L, an elec- 
tron unlikely keeps its phase coherence around a ring, 
therefore AAS oscillations are dominated by trajectories 
enlacing one ring only and all rings can be considered 
as independent. In the opposite regim^^H > tj^e 
interfering electronic trajectories explore regions much 
larger than the ring perimeter L. In this case, winding 



properties are more complicated (figure [T|) and the prob- 
ability Vn(t) may strongly differ from the one obtained 
for a single ring 'P™^{t). A theory must be developed 
to account for these topological effects, which leads to 
an harmonic conteiit qui te different from ([3|. This was 
done in Refs .1221251221211221 for large regular networks. This 
theory was later extended in Ref in order to deal with 
arbitrary networks, properly accounting for their connec- 
tions to contacts^. 




Figure 1; A square metallic network submitted to a magnetic 
field. Schematic picture of a closed diffusive trajectory wind- 
ing a flux 2(j) is represented ((f> is the flux per elementary pla- 
quette). 



(a) e-e interaction leads to geometry dependent 
decoherence. - Not only the winding probability Vn{t) in- 
volved in eq. ([2| is affected by the nontrivial topology of 
networks, but also (e'**-*-*) — e"*/"^^ describing the na- 
ture of phase coherence relaxation is replaced by a more 
complex function. When decoherence is due to e-e inter- 
action, the dominant phase-breaking mechanism at low 
temperature, this relaxation is not described by a sim- 
ple exponential anymore. This situation will be refered 
as model B throughout this paper. Such decoherence can 
be modeled in a one-electron picture by including dephas- 
ing due to the fluctuating electromagnetic fleld created 
by the other electron^. Therefore the pair of reversed 
interfering trajectories picks up an additional phase <i>[C] 
that depends on the electric potential V . Averaging over 
the fluctuations of the potential leads to the relaxation 
of phase coherence. The harmonics present the structure 



Act„ 



'v,cr , 



(4) 







where averaging is taken over the potential fluctuations 
{■■■)v and over the loops with winding n for time t, 
(• • •)cj"- In a quasi Id-wire, the relaxation of phase co- 
herence involves an important length scale, the Nyquist 
length Ln, characterizing the efficiency of the electron- 
electron interaction to destroy the phase coherence in 
the wire. We will see that, in networks also, the Nyquist 
length is the intrinsic length characterizing decoherence 
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due to e-e interaction. It is given bjEI^SISlISn] 



Ln — 



tfaoDs 



1/3 



1/3 



(5) 



where s is the cross-section of the wire. We have rewrit- 
ten the Nyquist length in terms of the thermal length 
Lt = \/hlT/ltBT, the elastic mean free path £e and the 
number of conducting channels (not including spin 
degeneracy) ; is a dimensionless constant depending 
on the dimension (ad = Vd/Vd-i where Vd is the vol- 
ume of the unit sphere in dimension dj^. In the follow- 
ing we will set h — ks ^ ^- In the infinite wire, the 
decaying function (e'*)y_Ct can only involve the unique 
length Ljy. A precise analysis of the magnetoconduc- 
tance (MC) of the wire shows that, in this case, the cal- 
culation of p| and Q for the infinite wire, for which 
V{t) cx l/y% leads to almost indistinguishable results 
provided \/2Lff. Therefore the analysis of the 
MC of the wire suggests that the sophisticated calcula- 
tion of Q can be replaced by the simpler one (pj) with 
Q-t/r^ ^ Q-t/2T^ ^Yiere tn = L%/D (x T'^/Ms the 



Nyquist time. However this is a priori not true any- 
more as soon as we consider networks with a nontrivial 
topology because electric potential fluctuations depend on 
the geometry, and therefore the decoherence is geometry- 
dependent. 

Let us formulate this idea more precisely. Being re- 
lated to the potential as ~ /* V, fluctuations of the 
phase <i> picked by the two reversed electronic trajecto- 
ries can be related to the power spectrum of the poten- 
tial, given by the fluctuation-dissipation theorem (FDT) : 



:x{t) where TZt 



J,,. )v,c, - e^TUt - ^x{t) where TZt = x{t)/{c7os) 
is the resistance of a wire of length x{t). The average 
(■ ■ ')v.Ct is taken over potential fluctuations and closed 
diffusive trajectories for a time scale t. The length 
x{t) is the typical length probed by electronic trajecto- 
ries. For an inflnite wire it scales like x{t) ^ \/ Dt there- 
fore {^'')v.,Ct - {t/rN?'^ where tjv = L'j^/D - T-2/3 
is the Nyquist time. On the other hand, in a ring, dif- 
fusion is constrained by the geometry : harmonics of the 
MC of a ring involve winding trajectories for which the 
length scale probed is therefore the perimeter x{t) ~ L, 
leading to (<i>^)y,Ct ^ (V'^'c) where Tc ~ T^^. Therefore, 
in a ring, depending on their winding, trajectories probe 
different length scales : Ljv oc T~^/'^ or Lc oc T"^/^. 

Let us summarize. At the level of eqs. (2[3l, 
is a phenomenological parameter put by hand. The 
modelization of decoherence due to e-e interaction of 
AAK shows that, in an inflnite wire, the WL correc- 
tion probes the Nyquist length Ljv oc T~^/^ (the only 
length scale of the problem). This shows that, in the 
MC of the inflnite wire, the phenomenological parame- 
ter must be replaced by L^^ — > (x T"^/'^. On 
the other hand the MC of a ring involves a new length 
scale Lc = L'^J'^ /L^^'^ combination of the Nyquist length 
and the perimeter. In this case, assuming the simple 



AAS behaviour A(t„ oc e "-^/^v^ the phenomenological 
parameter should be substituted by L^p —> Lc (x T~^/^. 

Geometry dependent decoherence in ballistic rings. - It is 
worth pointing that such a geometry dependent decoher- 
ence can also be observed in ballistic systems : potential 
fluctuations responsible for decoherence depend on the 
precise distribution of currents inside the device, that are 
affected by the way the current is injected through differ- 
ent contacts^. Depending whether the measurement is 
local or nonlocal, different phase coherence lengths have 
been extracted from the damping of AB oscillation^^. 
The different are probed by changing the contact con- 
flguration (current/ voltage probes )p^, whereas in the dif- 
fusive ring, the different length scales are probed by con- 
sidering different harmonics. 



III. GENERAL FORMALISM 

We flrst recall the basic formalism and apply precisely 
the ideas given in the introduction. We will consider the 
reduced conductivity a, deflned by 



2e' 



(6) 



where s is the cross-section of the wire. The reduced WL 
correction has the dimension of a length. As mentioned 
above, it is a sum of contributions of interfering closed re- 
versed electronic trajectories, which can be conveniently 
written as a path integral : 

Aa{x) = -2Pc{x,x) (7) 
= -2 dt P:j(r)e-^o'*"(3i^+2i-A(.))gi*[.(r)] ^ 

Jo Jx{0)=x 

Pc{x,x) is the so-called Cooperon. Summation over dif- 
fusive paths for time t involves the Wiener measure 
I?a::(T) exp — Jq dr (we have performed a change of 
variable t — > t/ D so that "time" has now the dimension 
of a squared length). Each loop receives a phase pro- 
portional to the magnetic flux 2 dr xA{x) intercepted 
by the reversed interfering trajectories, where A{x) is 
the vector potential. The factor 2 originates from the 
fact that the Cooperon measures interference between 
two closed electronic trajectories undergoing the same 
sequence of scattering events in a reversed order. Finally 
we have introduced an additional phase <i>[a;(r)] to ac- 
count for dephasing : dephasing due to penetration of 
the magnetic fleld in the wirea^ or decoherence due to 
electron-electron interaction. In this latter case the phase 
depends also on the environment dynamics over which 
one should average. 

The X dependence of Aa{x) in eq. In a general net- 
work, in the absence of translational invariance, the WL 
correction to the conductance was shown to be given by 
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an integration of Acr(a;) over the network, with some non- 
trivial weights attributed to the wires. Let us write the 
classical dimensionless conductance as g'^'^^^^ — adNcie/^ 
where the effective length £ is obtained from addition 
(Kirchhoff) laws of classical resistances (dimensionless 
parameter ad was defined above : 03 = 4/3, a2 = 7r/2 
and ai — 2). Then the WL correction to the conductance 



4301 



*^ . i J wire i 



(8) 



where the summation runs over all wires of the networks 
and li is the length of the wire i. Eq. ([8| was demon- 
strated in Rcf."^" for the conductance matrix elements 
of multiterminal networks with arbitrary topology. This 
result relies on a careful discussion of current conserva- 
tion (derivation of current conserving quantum correc- 
tions can be found in Refs.l^^ESllZIll^, This point will 
play a relatively minor role in the present paper. Eq. Q 
may be used in order to calculate geometry dependent 
pref actors. 

Magnetoconductance oscillations and winding 
properties.- In an array of metallic rings of same 
perimeter, the magnetic flux is an integer multiple of 
the flux (j) per ring dTxA{x) = 4> x Af[x{T)] where 
J\f[x{T)] g Z is the winding number of the closed trajec- 
tory (the number of fluxes encircled). This makes the 
WL correction a periodic function of the flux cf) = 9^, 
where 0o = h/e. The n-th harmonic of the MC 

r2^ d6l 



Act,, = 







27r 



Aa{e) e- 



(9) 



involves trajectories with winding number n. We can 
write the harmonics as 



Act„(x) EE -2Pi")(x,a;) = -2 / 

^0 



dt 



(10) 



/>x(t)—x 
lx{0)=x 



dr ^i<S>[x(T)] 



where the Kronecker symbol selects only trajectories for a 
given winding number n. Let us introduce the probability 
for a diffusive particle to return to its starting point after 
a time t, with the condition of winding n fluxes 



(*^5 ^) 



x{t) = 



x{0)=x 



For example, in an isolated ring of perimeter i, this prob- 
ability is simply given by 



'{x,x;t) 



Then, we can rewrite the harmonics as 



Adn{x) = -2 / dtVn{x,x;t){&^^^"^)c 



(12) 



(13) 



which is the structure given in eq. El). In eq. ( 13 1 we have 
introduced the notation C" = {xyr), r € [0,t \x{0) = 
x{t) — x ; A/'[x(t)] = n) for a closed diffusive path wind- 
ing n times. (• • designates averaging over all such 



paths, with the measure of the path integral (111. The 



phase <i>[C] accounts for dephasing and eliminates the 
contributions of diffusing trajectories at large time. We 
now discuss two possible modelizations for this function, 
denoted by "A" and "B" . 



A. Model A : Exponential relaxation 



The simplest choice is an exponential relaxation, with 
a dephasing rate 7 = l/r;p = 1 /L^ : 



(14) 



This simple prescription correctly describes d ephas ing 
due to spin-orbit coupling, magnetic impuritiea^^'^lL 
feet of penetration of the magnetic field in the wire a^, or 
decoherence due to electron-phonon scattering^^. Using 
( 12|13 1 with this exponential decay yields the familiar re- 
sult (|3j) for the isolated ring. 



Model B : geometry dependent decoherence 
from electron-electron interaction 



It turns out that the simple exponential relaxation 
does not describe correctly the decoherence due to 
electron-electron interaction, the physical reason be- 
ing that this decoherence is due to electromagnetic 
field fiuctuations that depend on the geometry of the 
sys tem. AAK have proposed a microscopic descrip- 
tiorP^H that we can rephrase as follows. In eq. Q, 
the phase <& picked up by the reversed trajectories de- 
pends on the environment (the potential V created 
by the other electrons due to electron-electron inter- 
action) : $y[Ct] = dr [y(x(r), r) - y(x(T), i - r)]. 
Averaging over the Gaussian fluctuations of V leads 
to (e'*^'['^'l)v = Q-ki'^vlCtf)v where the fluctuation- 
dissipation theorem (written for lo describing classi- 
cal fluctuations) {V{r,t)V{r',t'))v ^ ^S{t-t')Pd{r,r') 
gives l{<^>v[Ct?)v = T[Ct] t witlP 



T[Ct]t^-^ dTW{x{T),x{t^T)) 



Rk Jo 



dT'R{x{T),x{t-T)) , 



(15) 



(16) 



where Rk — h/e^ is the quantum of resistance. The 
function W{x, x') is related to the diffuson, solution of 

-APd{x,x')=5{x-x'), by 

^^^^ ^ Pd{x,x)+ Pd{x',x') _ p^^^^ ^^^^ 
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This function has a physical interpretation discussed in 
the appendix |E] : it is proportional to the equivalent re- 



sistance TZ{x, x') between the points x and x' (figure 21 1. 
With this remark, we see that eq. (16) can be under- 



stood as a local version of the Johnson-Nyquist theorem 
relating the potential fluctuations to the resistance. 

In eqs. we have introduced a decoherence rate 

r[Ct] which depends not only on the time but on the 
trajectory itself. Therefore the decay of phase coherence 
is now described by 



(e 



i^v [Ct 



IV fit 



-r[Ct]t> 



Ct 



(18) 



Within this framework, relaxation of phase coherence 
is not described by a simple exponential decay like in 
eq. (|14| but is controlled by a functional of the trajecto- 
ries^ x{t). Therefore the nature of decoherence depends 
on the network, through the resistance 'R-{x, x') between 
X and x' , and on the winding properties of the trajecto- 
ries. 

The central problem of the present paper is to compute 
the path integral 

/■OO j*x{t)—X 

Aan{x) = -2 die-T* / Vx{t) 5,t.Af[x(r)] 

Jo Jx(0)=x 

-J^dT [lx^+-^W{x{T),x{t--r))] 



X e 



(19) 



for the different networks. Such a calculation has been 
already perfor med in two cases : the infinite wir^ and 
the isolated rin^^. 

The logic of the following sections is the following : 
first we study the winding properties in the network. 
For that purpose we first compute the WL correction 



within model A, eq. ^l^ with (flil. The probabil- and 



ity 'Pn{x,x;t) can be extracted from an inverse Laplace 
transform with respect to the parameter 7. Having fully 
characterized the winding properties, we use this infor- 
mation in order to study the harmonics Act„ within the 
model B describing decoherence due to electron-electron 
interaction, eq. (13 1 with (18 1. 



IV. THE WIRE AND THE RING 

We first recall known results within the framework of 
model B concerning the simplest geometries that will be 
useful for the following. 



the path integral 



(20) 



-2 / die-"^* 

. .-So<^r[\x- + ^\xir)-x(t-r)\ 

Jx{0)=x 



can be computed thanks to translational invariance (as 
pointed in Ref.^-'^, using the symmetry of the path integral 
we can perform the substitution x{t) — x{t — t) x{t), 
provided that the starting point of the path integral is 
set to X —> 0, see appendix [A| . Combining exponential 
relaxation (model A) and decoherence due to e-e interac- 
tion {model B) allows to extract the function ( 18 1 with 
an inverse Laplace transform of the AAK resullP'^lSlIll] 



A£^: 



dt 



= L 



N 



yi^t 
Ai'(7i^) ' 



(21) 
(22) 



where Ai(z) is the Airy functiort ^. As mentioned 
above, this expression is very close tcP^H Act ~ — ( -I- 



1 

2Ll 



-1/2 



the result obtained by performing the substi- 



tution (e'*^['^'l)y,Ct e-*/2i'jv (see figure |2|. 

The inverse Laplace transform of ( pT| was computed 
in Ref.'^ with residue's theorem : 



/wh-c(i/Tw) = (e'*^'['^*l)yA (23) 



m— 1 



(24) 



where Um are zeros of Ai'(z). In particular ux ~ —1.019 



(m — |)]^/^ for m 00. The limiting 



behaviours are 

/wire (t/TN) — 1 



/TT / t 

4 \tn 
\ui\ V Tat 



3/2 

for t < TN (25) 
Q-\u,\t/r^ forrAr<t. (26) 



Note that the short time behaviour can be obtained by 
expanding / (t/rAr) - (e'*)y,c. = (e-^C^')^)^^ ~ i „ 
\{'^v[Ct\^)vfif This limit can be simply obtained by 
noticing that in the wire W{x{T),x{t — r)) ~ x{t) ^ \/t, 
therefore V\Ct\t - Tt^l"^ - (tjTf^fl'^ where we recover 
that the Nyquist time scales as rjv oc T~^l^ . 



Phase coherence relaxation in an infinite wire 



Ref.[3]. In this case we have Wwiie(a;, a;') 



\x — X 



B. Phase coherence relaxation in the isolated ring 



For the isolated 



The case of an infinite wire was originally solved in Wring(a;, x') 



ring 



of perimeter L, we have 
— i?— i^). The path integral 



and (19 1 can be computed exactlj^ (see appendix [a|) . Up to 



LIU 



2 3 



Figure 2: Comparison between —hi{x^)/Pd'{x'^) (black con- 
tinuous line) and 1/ \/l/2 + x'^ (red dashed line). Relative 
difference (inset) does not exceed 4%- 



a dimensionless prefactor, we obtain 

i(L/L«)='/" 



for Ln <.L (27) 
(28) 



This result can be simply understood as follows : in 
the ring, trajectories with finite winding necessarily ex- 
plore the whole ring. This "ergodicity" implies that 
W{x{T),x{t — t)) ~ x{t) ~ L and therefore the deco- 
herence rate V\Ct\ t ~ TLt ~ t/rc involves the different 
time scale Tc ^ 1 / {TL) , according to the physical ar- 
gument given in section [ll] As a consequence Eq. ( 13 ) 
indeed leads to- 



where 




(29) 



(30) 



Figure 3: The relative difference (1 — ri{x^)/ ^^{Tv/Sy + x^) 
between the effective length L'^'^ , eq. (33 1, and its approxima- 
tion, eq. (361, does not exceed 1.5%. 



AtT„ oc fn{L/L^). Eq. (27 1 gives fn{x) — G 



,3/2 



with 



the geometry independent phase coherence length L^p 



The temperature dependence Aa„ 



_ ^3/2^1/2 



was 



first predicted in Ref.El using instanton method (with a 
different pre-exponential dependence) and studied in de- 
tails in Ref.'^ where the path integral (19 1 was computed 
exactly for the isolated ring. The effect of the connect- 
ing arms was clarified in Ref.'^. The fact that the pre- 
exponential factor is is related to the fact that the 
smooth part of the MC, due to the penetration of the 
field in the wire, probes the same length scale as in the 
infinite wire. 

It is worth pointing the recent work'^ in which the 
crossover to the Od limit is studied in a ring weakly con- 
nected. In this case the authors get a crossover from 
cx (diffusive ring) to oc (ergodic) for 

temperature below the Thouless energy. This latter be- 
haviour coincides with the result known for quantum dots 
in the same regimd^. 



Phase coherence length : L^p oc T or L^p oc 
7^—1/2 ?_ jvj^ote that the introduction of a new length 
scal^Lc might appear arbitrary since the harmonics may 
be written uniquely in termd^of L and Ljy. The differ- 
ence between ([27| and ( 29 1 is a matter of convention and 



may be related to the experimental procedure. The usual 
method extracts the phase coherence length from the 
analysis of MC harmonics. Then it is natural to see how 
the winding number n scales with the phase coherence 
length, or more properly how the length nL scales with 
and therefore assume the form Aijn oc f{nL/L^). 



From eq. ( 29 1 we see that the function is simply the ex- 
ponential, /(i) — e~^, with a perimeter dependent phase 
coherence length ^ Lc (x (TL)~^/^. Another proce- 
dure may consist in studying the harmonics content as 
a function of the perimeter L, that is for different sam- 
ples. The experiment is then analyzed with the form 



C. Penetration of the magnetic field in the wires of 
the ring 

Networks are made of wires of finite width w. The pen- 
etration of the magnetic field in the wires is responsible 
for fluctuations of the magnetic flux enclosed by trajec- 
tories with the same winding number but different areas. 
In the weak magnetic field limit, this effect is described 
by introducing an effective dephasing rat^^^l 



eBw 



(31) 



The question of how to combine the two decoherence 
mechanisms {models A & B) in the ring was discussed 
m Ref.'^. It was shown that the WL correction of the 
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ring presents the structure 

Ao-„ ~ Ln 
for Lf4 <^ L, with 



Ai'(L^/L|) 



(32) 



(33) 



where ri{x) = (j 



x) arctan ^ 



2^ 



Pref actor. - In eq. ( 32 1 , the pre-exponential factor co- 
incides with the result obtained for an infinite wi rj^. 
The ratio of Airy functions can be approximated a^^Ill] 
-^^r^ — — (^ +a;)^^/^ (figure 2j). In other terms, we may 
write the zero harmonic (i.e. the result for the infinite 
wire) as 



03 




Figure 4: Typical shape of the MC curve of a network (here 
a cham of rings). Rapid oscillations are AAS oscillations. 
Damping of oscillations over large scale is due to the penetra- 
tion of the magnetic field in the wires. 



where 



T cnv 



21/^ 



r env 



r2 



-1/2 



(34) 



(35) 



This combination expresses that, in a wire of width w, the 
penetration of the magnetic field provides the dominant 

cutoff when typical trajectories enclose more than one the MC curve. Which {Ln or L^) is obtained from 
quantum flux L^wB > (f>o (here - for trajectories g^^h a curve ? According to ^ we see that -> Ln 
with winding n — 0). 



curve of a network. The following discussion applies 
to the case < L where (32 1 holds. Figure |4] rep- 
resents a typical MC curve, here for a chain of rings. 
It exhibits rapid AAS oscillations with a period given 
by ;Bosc ^ 4>o/L'^i superimposed with a smooth varia- 
tion over a scale ;Bdamp ~ (p^/wLip. The phase coher- 
ence length can be extracted either from the amplitude 
of the oscillations or from the decay of the envelope of 
the MC curve. Which {L 



for 



Exponential damping.- In the exponential of eq. (32 1 
the effective length interpolates between L°^^ — §i, 
Lb > Lc and L"^" ~ Lb for Lb < i, 
haviour is well approximated by 



osc 

Its overall be 



r osc 



8L, 



1 



-1/2 



1 

^B 



-1/2 



(36) 



which differs with (33 1 by less than 1.5% (figure[3|. When 
Lb is the shortest length, the decay of AAS oscillations 
~ g-rai/i-B (2an be understood from the fact that modula- 
tions of the fiux enclosed by trajectories with finite wind- 
ing become larger than the quantum flux nLwB > 0o. 
Eq. ^36] ) was used in the analysis of the recent experi- 
meniPr 

These two remarks show that the magnetic length Lg 
probes two different length scales : in the pre-exponential 
factor Lb probes the Nyquist length Ln oc T^-^^^, 
whereas in the ratio of harmonics Acr„/Aao, the mag- 
netic length Lb probes the length scale L^ oc T^^/^. 



D. How to analyze MC experiments in networks 

In order to understand the implications of this re- 
mark, let us discuss the structure of the typical MC 



Nn '^■^ AAK theory (T"'"^) 
.-0.34- x2 



T 



1 »«0 s 



o <p 



-r-0.53 Q, 



0.1 



1 



T(K) 



Figure 5: Phase coherence length as a function of the temper- 
ature obtained from measurements realized on a large square 
network of lattice spacing a — 1 /im (perimeter L — 4 fj,m). 
The length scale — > LSf" oc T"^/^ has been extracted from 
the smooth MC while the length L^p L"^"^ oc T~^^^ has been 
extracted from ratio of h armo nics . (Note that L™^ and L'^'^ 
in Ref^ denote eqs. (35 I and (36 1 for Lb = oo). 



in the pre-exponential factor, which mostly dominates 
the smooth envelope, while L^ —^ L^ in the exponential 
decay, which dominates the damping of the rapid oscil- 



lations. In order to decouple the two effects the analysis 
of t he ex periments of RefP have been analyzed as fol- 
Jq^Ji6|65| . |.]^g Fourier transform of the MC curve Aa{B) 
presents broadened Fourier peaks due to the penetra- 
tion of the magnetic field in the wires. Integration of 
Fourier peaks eliminates this effect. Ratio of harmonics 
involve the length scale Lc- The length L^v was extracted 
from the smooth envelope AcTenv = J2n '^^n — Actq for 
Ln L. The temperature dependence of the phase 
coherence length was extracted in this way in Ref.l^. Re- 
sults are plotted on figure [5] exhibiting clearly the two 
length scales in the regime Ln ^ L. We see that it is 
crucial to analyze the experiment in terms of the MC 
harmonics Acr„. 

Isolated ring vs ring embedded in a network. In 

transport experiments the ring is never isolated : it is 
at least connected to contacts through which current is 
injected. Moreover the samples are often made of a large 
number of loops, in order to realize disorder averaging. 
The results obtained for the isolated ring are fortunately 
relevant to describe a more complex network of equiva- 
lent rings (Figs, [l] & [9]) when the rings can be consid- 
ered as independent, i.e. when interference phenomena 
do not involve several rings ; this occurs when <^ L 
(or Ln <C L), in practice in a high temperature regime. 
This temperature dependence of harmonics is rather dif- 
ficult to extract from measurements since harmonics are 
suppressed exponentially. This has been done only very 
recently in Ref.l^. Another difficulty is that the "high 
temperature regime" Lpf L is in practice quite nar- 
row in these samples due to fact that electron-phonon 
interaction dominates the decoherence above 1 K (in the 
sample of Refs.l^^^^Ljv is much larger and when <^ L 
the role of electron-eletron interaction is negligible). 

It is an important issue to obtain the expression of the 
WL correction for a broader temperature range, that is 
to study the regime L < Lj y. This regime is reached 
in several experiment d^^ * ^^ *^. In this case diffusive inter- 
fering trajectories responsible for AAS harmonics are not 
constrained to remain inside a unique ring, but explore 
the surrounding network (see figs, [l] [6]& |9|. This affects 
both the winding properties and the nature of decoher- 
ence. The MC oscillations are therefore network depen- 
dent. In the following sections we discuss the behaviour 
of the MC harmonics in the limit L < (or L < Ljv) 
for different networks : a ring connected to long arms, a 
necklace of rings and a large square network. The case 
of a long hollow cylinder will also be discussed. 



THE CONNECTED RING 




Figure 6: A ring connected to reservoirs through long wires a 
and b. In the regime L^p > L, the WL correction is dominated 
by trajectories exploring the arms, as represented here. 



Model A 



Let us consider first the case <^ L. The Cooperon 



is constructed in appendix F 1 (see Ref.l^ and we obtain 



Aanix e ring) 



-nL / 



(37) 



where x is any position inside the ring far from the ver- 
tices (at distance larger than L^p). 

• Now we turn to the regime L^p ^ L. The Cooperon is 
uniform inside the ring and is computed in appendix |F 1[ 
We have^ : 



for L < (38) 



where x is any position inside the ring, or in the arms at 
a distance to the ring smaller than L^p. We emphasize 
that this behaviour, quite different from (|3|, is due to 
the fact that the diffusive trajectories spend most of the 
time in the wire^ (the distribution of the time spent by 
winding trajectories in the arm was analyzed in Ref .'^ . 

Winding probability in a ring connected to 
arms.— We now derive the winding probability for a ring 
connected to Na infinitely long arms (figure ItI) from the 

inverse Laplace transform of the Cooperon Pc:'^\x,x) — 
— ^Adnix). At small time, (37l gives 



T^ji (x, X , t) 



Na 



1 



for t < (39) 



where x is inside the ring, far from a vertex (at distance 
larger than 

At large time scales, t ^ L^, the arms strongly mod- 
ify the winding properties around the ring : the time 
dependence of the typical winding number becomes sub- 
diffusive nt ~ t^/^, to be compared with the behaviour 
for the isolated ring rit ~ t^^'^ reflected by eq. (12 1. For 



a ring connected to Na inflnite wires, eq. (12) is replaced 
by the probability^ 



In this section we consider the case of a single ring 
connected to two wires supposed much longer than Lp 
(flgure[6|. This problem has already been considered in 
Refs.5 ^. 
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where Na is the number of arms (as far as x is in- 
side the ring or at a distance to the ring smaller than 
^1/2 l/^ = L^, the Cooperon, or the corresponding 
probability is almost independent on x). The function 
^-(0, given by7 61 



= - Re 



du^/ue 



(41) 



is studied in the appendix[B]and plotted in the conclusion 
(figure 17). 




Figure 7: A ring connected to Na 



The general expression describing the crossover between 
(43 1 and (44 1 can be obtained easily using the formalism 
of Ref.30. 



B. Model B 

We now compute the harmonics of the conductivity 
within model B. We have now to consider eqs. (|13|15|18|). 
The function W{Xtx') has been constructed in Ref.^. In 
the limit Z^, li, ^ L and if x and x' belong to the connect- 
ing wires for x, x'(i.e. iAr)<C la, h, the function coincides 
with the one of the infinite wire W{x,x') ~ ^\x — x'\. 
Therefore, since in the limit L Ljq the diffusive tra- 
jectories spend most of the time in the wirePi' (figure |6|, 
the dephasing mostly occurs in the wires and the re- 
laxation of the phase coherence is similar to the one 
for the wire, eq. (24 1, irrespectively of the winding 



variable |u, 
the conductivity 



ml)- 



Introducing (24 1 in (13 1 and performing the change of 

\t/TN = 



we obtain for the correction to 



From the conductivity to the conductance. In 

the geometry of figure [6] the conductance is not simply 
related to the conductivity. The classical conductance 
of the connected ring is given hy g = a^lSS ^.(-el with 
L, — la + lc\\d + 4 where li is the length of the wire i. 

is the equivalent length. From eq. 



c\\d c 



^9n = 



{la+lc\\d + lb)'^ 



a {Ic + ldYJc {Ic + ldYJd Jb 



(42) 

dx Act„ (x) . 



The Cooperon Act„(x) = —2Pc {x,x) has been con- 
structed for different positions of the coordinate x in 
appendix Fl Depending on the ratio L^/L, the WL 



correction A(7„ is dominated by different terms. 



For L„ <C L, eq. (F6l shows that harmonics of the 



Cooperon decay exponentially in the arms (figure 24 1 ; 
inside the ring, the Cooperon ( |F8[ ) is almost uniform, 
apart for small variations near the nodes. Therefore A(7„ 
is dominated by integrals and in the ring and we 
have 



A5„ 



,2-Aa-„(x G ring) 



for L , 



(43) 



where x is any position inside the ring far from the ver- 
tices (at distance larger t han L^ ,). 
• For Lip ^ L, using eqs. ( F6|F8 1 we see that the terms 
and bring a contribution proportional to the perimeter 
L whereas the terms and bring larger contributions 



proportional to 
therefore 



J^dxP^"\x,x) 



iL^Pi")(0,0), 



-^Adnix e ring) for L <^ 



(44) 



Aa„(x) ~ - V2nLNL ^ 

x^"d«v^yvi/|^|u^|i/4„y^ 



(45) 



L 



N 



where x is any position inside the ring. Using (41 1 we 



rewrite the double integral in polar coordinates and per- 
form integration over the radial coordinate. We find 



Aan{x) ~ —Ln Fi (n^/2L/LN 



(46) 



with 



2^/2 

m—1 



and 

5(A) = Re 
1 

"71 



Re 



(48) 
(49) 



A convenient representation can be obtained by a rota- 
tion of 7r/4 of the axis of integration in the complex plane. 
We get 



-A r°° 



2 o-At 



(50) 



The function i^i(^) is plotted on figurejsj We now analyze 
the limiting behaviours. 
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Figure 8: Acr„ oc Fi{^) (continuous blue Itne) as a function 
o/l/C^ oc Lf^/n^L. Comparison with the limiting behaviour 
<C 1) ~ 1.191/,f- 0.877 (red dotted line). Inset : ^f'i(C) 
in semi-log scale. For moderate values of ^, the function is 
well fitted by 1.21 e""'^* (magenta dots). Green dashed curve 
is the asymptotic expression, eg. (51 1. 



We first consider the limit Ln <^ n L. Eq. (50 1 gives 



^(A»i)c.-l5e-Aj»_^(l_e-A.) 



4\/2 



VttA 



therefore Fi{i > 1) ~ i|ui|-i3/8y^e-l"il « and 



function but only to the scale over which it decays. In the 
other limit Lat <C n^L, the calculation of the harmonics 
rather involves the tail of the function (e'*)^^^. Eq. (51 1 
presents an exponential decays similar to eq. ( 38 1 , with 
a different pre-exponential power law since the decay 
(e'*)y,Ct « is different from the simple 

exponential decay e~*/'^»' for model A. The additional ^/t 
in model B explains the different pre-exponential terms 



in eqs. (38) and (51 1. 

Could we map the results of the two models through a 
simple substitution of phase coherence length, as for the 
infinite wire ? In the regime L <C <C ri^L we should 



compare the exponentials of (38 1 and (51) what leads 
LnI\/W 



to L,, 



I ~ O.OQLat, however pre-exponential 
factors cannot be matched, obviously. In the regime 
Lm ^ n^L we rather compare the square roots (38 1 and 
(53) and therefore A\Lm — 1.418LAr. Despite 



there is no unique simple substitution, we get in both 



cases L. 



Conductivity vs conductance. We discuss the rela- 
tion to the conductance. In the regime Ln ^ L discussed 
in this section, we expect 



A5„ 



£2 



Act„(x) 



(54) 



for X G ring, for the same reason as for model A, eq. (44| 
(the factor comes from the fact that the integral (Isl 
is dominated by regions of typical size L^r in the arms) 



VI. THE CHAIN OF RINGS 



Ao-„(x) 



'2t|Ml|T 



1/4 

rfiL 



(51) 



_ y-i/i2g-„ lV^tV« fori«L^«„2^, (52) 
where = y/2\ui\^/^ ~ 1.421. 



For the lowest temperature n^L ^ Ljv, eq. (49) gives 
5(0) = ^r(3/4)2. Therefore Fi{^ « 1) ~ 1^/^ with 

Ai = i^/i7^^(3/4)2^^^l \urn\-'/^ ^ 1.191 (the sum 
is ^ 1.989) 



Aan{x) ~ -Ai 



LnL 



for n^L < L 



N 



(53) 



(a) 




Figure 9: Chains of rings. If we consider the regime la ^ 
L^{or Ln) 3> L the rings of figure (a) can be considered as 
independent. On the other hand, in the regime L^p{or L^) 3> 
la, L interferences occur between trajectories encircling sev- 
eral rings ; these latter cannot be anymore considered as in- 
dependent in this case, as illustrated by the trajectory m the 
network (b). 



Comparison between models A &; B. We have seen 
that for the wire, the MC obtained from the t wo m odels 
are related through L^p \/2 Ln {cf. section IV I . It is 



tempting to look for a similar relation for the connected 
ring in the limit L <C Lj^. 

Let us compare the results for the two models of de- 



coherence. In the limit n L ^ L^f the expressions (53) 



is very close to (38) because in this case, the harmonics 
involve an integral over time of the function (e'*)y.Cf 
Therefore harmonics are insensitive to the details of this 



Let us consider now a chain of rings separated by arms 
of lengths la (figure [9] a). The case where the phase co- 
herence length remains smaller than the arms (L^ < la) 
can be obtained from the results of the previous section 
since the rings can be considered as independent. The 
conductance of a chain of such iV^ rings is given by per- 
forming the substitution , nr , — ra- in the 

expression of Ag„ for one ring, eqs. ( 43|44|54 1. 



When arms separating the rings are smaller than the 
phase coherence length {la < L^), coherent trajectories 
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enclose magnetic fluxes in several rings. In order to study 
how the MC harmonics are affected by this effect we con- 
sider the limit when rings are directly attached to each 
other (figure [9] 



A. Model A 

Considering the conductance of the symmetric chain 
of figure [9]b, the weights of the wires involved in 



eq. ([8| are all equal. This justifies a uniform integra- 
tion of the Cooperon in the chain. In this case we 
can use the relation between the WL correction and 
the spectral determinanlP^l2ll2ll (appendix [pj). The spec- 
tral determinant of the infinite chain is given in ap- 
pendix |F 2[ Averaging the Cooperon in the chain, 
Aa{0) Er7^^.^^Aa(x,0), and using ( |D5|F10| , we fi- 
nally obtain^ 



' coth(L/2L,) - ^ + ^ sinh(X/2L,) 

^ ^cosh^(L/2L^) - cos2(6'/2) / 



(55) 



where 9 = ^^(p/cpQ is the reduced flux per ring. We now study the harmonics rt ^ : 



2ir 



d0 sinh(i/2L^)e'^ 



V2 Jo 27r ^cosh(L/L^) - cos(6i) 



-L^ smh{L/2L^) 



dz 



cio 2i7r ^(e^/i^ - z){z - e-^/^-p)/z 



(56) 



The integration in the complex plane is performed along the unit circle in the clockwise direction. The segment of 
the real axis [0,e~^/'^'^] is a branch cut. The contour of integration is deformed to follow closely this segment. We 
obtain : 



Aa„ 



L sinh(L/2Ly) ^_^„-i/2)l/l^ 
2tt L/2L^ 



1 , (l-u)"-i/2 
au 



\/v? + eu 

with e = q2L/l,^ _ -y^g recognize the integral representation (CI) of the hypergeometric function*^ F 



(57) 



AcTn 



L sinh(L/2L^) 



27r L/2L^ 



fl g-(n+l/2)L/L^, ^ 



- ,n-\ — 
2' 2 



F 



- ,n-\ — \n 
2' 2' 



(58) 



where B(x,y) is the Euler /3 function. 

Weakly coherent limit.- We consider the limit L^p <C L. 
Using F{a, /?; 7; ^ ^ 0) ^ 1, we obtain : 



(2n- 1)!! ^ 



iL/L^ 



for L^<t: L, (59) 



a result reminiscent of the result of the isolated ring ([s]) , 
with a different prefactor originating from the probability 
to cross the vertices of coordination number 4 (note that 

i^n^m _ r(n+l/2) 1_ r „ ^ 1 ^1 

Large coherence length. - In the opposite limit limit ^ 
L. Eq. (|C4| gives : 



for L-€.L,^<^nL. (60) 



We have recovered an exponential damping of the har- 
monics, reminiscent of ( 3|59 ), but with a different 
dependence of the pre-exponential factor. 



On the other hand, for harmonics with nL <^ L, 



the harmonics can be expanded by using eq. ( |C3| . Let 
us introduce 6„ — In n — ^(n + ^) -t- V'(l)- These co- 
efficients converge to a finite limit at large n : bn = 
-C-24^ -t-0(;ij-), where C= -^'(1) = 0.577215... is the 
Euler constant. Finally we obtain 



ABr. 



L 
'2n 



[ln{2L^/nL) + 6„] for nL ^ . (61) 



It is useful to remark that the expressions (60 1 and 



(61 1 coincide with the limiting behaviours of the modified 



Bessel function Kq^z) for large n (the proof is given in 
appendix O : 



Act, 



L 

'2^ 



Ko{nL/L^) for L <Si. L^ 



(62) 



Up to a factor 1/2 interpreted below, this expression coin- 
cide s wit h the MC harmonics for a long hollow cylinder'^, 
eq. (102) recalled in section VIII We compare this ap- 



proximation with the exact expression ( 58 1 on figure [10 
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We see that the approximation is already good for n = 1 , 
provided L^p > L. The difference rapidly diminishes as n 
increases. 



b 
I 




Figure 10: Comparison between exact result (581 for n 



(dashed blue line) and 2 (dotted red line) and the approxi- 
mation (621 (black conttnuous line). The interrupted green 
curve is ( |59| |. Even for small n (= 2) eq. ( |62[ ) is a very good 
approximation of eq. ( 58 1 . 



that this probability is similar to the one found for an 
infinitely long hollow cylinder, apart for a factor 1/2. 
This additional factor can be understood from the fact 
that, starting from a given ring, it is equiprobable to 
return in one of its two arms. 

Let us give a heuristic argument to recover roughly 
(64 1, that will be useful for the following. Arriving at 



a vertex, the diffusive particle equiprobably chooses one 
of the four arms. Therefore it is equiprobable to wind 
a ring or not, while diffusing along the chain. This sug- 
gests that the winding probability is almost independent 
on n, up to n '--^ Vi/L, the maximum number of rings 
explored for a time t. This rough approximation would 
be Vn{x,x;t) ^ NtO{\/t/ L — \n\). The normahsation is 
estimated easily : since diffusion along the chain is one- 
dimensional, we expect 7'„(x, a;; ~ 1/v^ so that 
Vn{x,x;t) ^ L/t for \n\ < \JtlL and otherwise. This 



is a crude estimate of eq. ( 64 1 



B. Model B 



Logarithmic divergence of the harmonics for 



oo. 



We see from eq. (61 1 that the harmonics are 
weakly dependent on n (for n <C L^pjV). This logarith- 
mic behaviour reflects the singular behaviour A(Tosc(^) ~ 
1/61, with a cutoff at 6* - LjL^ : Aa^ J^^l^ f = 
ln(L(p/(nL)). The harmonics are therefore almost inde- 
pendent on n as soon as n is small enough compared to 
Lip/L. Note that in practice, this logarithmic divergence 
of the harmonics is limited : when the phase coherence 
length reaches the total length of the chain L^^, ~ N^L, 
harmonics reach a finite limit, AUn ~ — ^ ln{Nr/n), due 
to the effect of boundaries (external contacts) . 

Winding probability.— We now extract the probability 
Vn{x,x\t) from these results. First of all the behaviour 
( 59 1 is related to 



T^n{x, X, t) 



(2n-l)!! 1 
2"+in! VIt 



for t < 



^ (63) 

We have recovered (12 1 with an additional dimensionless 



factor coming from the probability to cross the vertices 
of coordination number 4 (this factor can be understood 
when one writes the trace formula for the heat kernel in 
the networkp2168|). 

The regime > L for the WL correction probes the 
regime t~> for the winding probability. We use the ap- 
proximation ( 62 ) in order to perform the inverse Laplace 



transform. Using the integral representation of the mod- 
ified Bessel functiorP^^ we get : 



T^n (*^; -^5 ^) 



L 

87rt 



for t > 



(64) 



We may check that ( 60|61 1 coincide with the limiting 
behaviours of this probability. It is interesting to point 



In order to compute the MC harmonics we first need 
to construct the function W entering in eq. (151. Fol- 
lowing appendix |f] we introduce a coordinate x — {x, /) 
to locate a point in the chain (the continuous variable x 
measures the distance along the chain while the discrete 
index / g {m, d} precise the arm, up or down). If x and 
x' do not belong to the same ring we have 



W{x,x')^^\x 



_ 1 

— 2 ^^wirc 



{x,x') 



(65) 



Remembering that W{x, x') is proportional to the resis- 
tance between points x and x' this equation has a clear 
meaning : when two consecutive nodes are linked by two 
wires instead of one, the resistance is diminished by a fac 
tor of 2 
to 



In the limit t :s> Lj^ the trajectories contributing 



'V,C" 



= e 



- J^dTW{xir)Mt-T)) 



(66) 



are trajectories extending over distances Ln ^ L along 
the chain. In this case we can neglect the contributions 
to the integral where the two arguments of W{x,x') are 
in the same ring. We have seen that for nL < \ft the 
measure of the Brownian paths weakly depends on n, 
therefore we expect that (e'*)y_Ct" — (G'*)v,Cti where the 
average of the l.h.s. is realized among Brownian curves 
of definite winding whereas the average of the r.h.s. is 
among all Brownian curves. The argument shows that 
the function describing decoherence corresponds to the 
result of the infinite wire in which L% 2L%. This 
factor 2 stands from the ratio between the resistance of a 
wire and of a chain of rings of the same length. Finally : 



(e'*-[^"V,Cr ^ /wirc(t/22/V^) 



(67) 
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We can use (24) in order to compute Ao-„ ~ 

-2/l°? dt UUt/2'/'L%). The lower cutoff 

takes into account the fact that expression of Vnix, x; t) 
is only valid for t ^ L^. However, except if n = 0, the 
cutoff is not important and can be replaced by 0. Using 
( 24 1 we obtain : 



-LF, 



where 



1 °° 1 



4 |Wr, 
m— 1 ' 



|3/2 



(68) 



(69) 



In the limit L jv <C nL the first term of the series dom- 
inates : 



Ao-„ ~ - 



L 



4|lii|3/2 



- K^nh j h ^ 



-nLT 



1/3 



e""'"" for L < La, < riL , (71) 
where K3 = 2-i/3|y^|i/2 ^ q.SOI. 




2 3 4 5 6 7 

1/|(i.e. VnL) 



Figure 11: A(t„ oc -^2(0 '^•5 <^ function 0/ 1/^ oc Ln /nL (blue 
continuous line). Comparison with the limiting behaviour (red 
dotted line) ^[ln(l/0 + Ccy\] with Ccyi ~ 0.51. 

In the opposite limit, for nL <C L^v (^ ^ 0), we can 
replace the sum by an integral and use the asymptotic 
behaviour Um ~ — [^|^(m — |)]^/^ 



5^0 bTT m 



(72) 



1 /"^/^^ dm 1 

- / — + cste = - [ln(l/e) + C^yx] . (73) 

DTT J I m ZTT 

Finally we find a result similar to the one obtained for 



exponential relaxation of phase coherence : 



Ao-, 



L 



[ln(Ljv/nL) + Cchain] 



for nL <C L 



N ■ 

(74) 

The constant is estimated numerically : we find Ccyi ~ 
0.51, hence Cchain — 0.74. This result could also 
have been more simply obtained by noticing that 



/wirc(i/2^''^L^) cuts the tail of 7'„(x, a;; t) on a scale L 



N 



(Ao-,, 



-2 f, 

J in 



(nL) 



. dt 



S-TTf 



Comparison between models A &: B. As we have 
done for the infinite wire and the connected ring, we es- 
tablish some correspondence between the results for the 
two models when Ljv S> L. 

In the regime L ^ Ljv ^ n^L the exponentials of (60 1 
and (70) may be matched if L,^ — > Ljv/'^s — 1.25L 



7 AT. 



In the regime Ljv 3> n L the logarithmic behaviours 



(70) (jei) and ([741 coincide for L,^ -> 1.87 L 



N- 



For not too large n, the two curves Act„ (Ln) w 
1.87iAr) are very close, apart for Ln <C L 



Aai^\L 



for which there is a qualitative difference between ([s]) 
and (pTl). 



From the conductivity to the conductance.— The 

weights in eq. ([S]) are all equals and conductance is re- 
lated to an uniform integration of Aa{x) in the chain 
The dimensionless conductance is given by Ag 
where Nr is the number of rings of the chain. 



VII. THE SQUARE NETWORK 

The easiest way to realize disorder averaging experi- 
mentally is to use netw orks with a lar ge number of rings, 
like 2d ne tworks (squarePnUIIIlIlIIll^ honeycomlPEaMl^ 
(j^g (jii | 5i |^ The "high temperature" regime {Ln L) is 
now well understood theoretically and experimentalljEl, 
but low temperature experimental results are still unex- 
plained'^^. Therefore understanding the magnetoconduc- 
tance of large networks when dccohcrence is dominated 
by e-e interaction still deserves some clarification. In this 
section we study the case of an infinite square network of 
lattice spacing a (figure [T]) . 



A. Model A 

The weak localization correction was derived analyt- 
ically by DouQot & Rammal (DR) for rational fluxes 
: 2TTp/q with p, q £ N (reduced flux is defined as 
2(j)/(po). They obtained^ 



27r 
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Aa(6'p,g) = - 



coth 



a 



i sinh - 

nq 



P;^(4cosh£ ) 
Pp,g(4cosh 



K 



Pp,g(4coshj^) 



(75) 



where Pp^q{ e) is a polynomial of degree q defined in ap- 
pendix G 3 where derivation of (|75|_2S recalled. K(a;) is 



the elliptic integral of the first kind^. 



Weakly coherent network. - The harmonics are suppressed 
exponentially as Act„ (x —LipQ~^^°'^^'f. Despite there 
is no close expression of the remaining dimensionless n- 
dependent factor, a systematic expansion of the spectral 
determinant can be written thanks to the trace formula 
of Ref)^ (the first terms of this expansion are available 
in RefP). 

Large coherence length. -The rest of the section is devoted 
to the large coherence length regime L^p ^ a. 

Continuum limit.— In the limit of small flux, <^ \, 
and large coherence length, ^ a, the discrete charac- 
ter of the network disappears and one should recover the 
results for the 2d plane in a uniform magnetic field. In- 
formations can be extracted from the study of this limit. 

The zero field WL correction is obtained from eq. (75) 
with p = q = 1, using Pi J ( x) = — x. Us ing the expansion 
of the elliptic integraP^SlI we fincP^ol . 



Aa(0) ~ ln(4L^/a) + 



6 



(76) 



This result is reminiscent of the WL correction of the 
film (91 1, but here, the cutoff at small scales is naturally 



provided by the lattice spacing a. 

The limit of small fluxes is studied in details in ap- 
pen dix IG 2[ Using that = 47r(/i/0o = 4:TrBa'^ /(j)o, 
eq. (|G4| reads 



A5-(6i< I)-Act(O) 



a 

2^ 



Am 



(77) 

This expression gives a quadratic behaviour for small 
flu^P 



Act(6') - Act(O) 



for 6* < (78) 



487r \ J 

and a logarithmic behaviour for intermediate fluxe^^ 

2 

Aa(0) ^^[lne + Csn] for « « 1 , (79) 
271- 

where Cs„ = -C - 3 In 2 - f ~ -3.704. 

We now turn to the analysis of the MC harmon- 
ics. A first simple remark allows to get the scaling 
of harmonics with time : the reduced fiux 9 is the 
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Figure 12: WL correction for L^/a — 2 (blue diamonds), 
L^p/a — 5 (brown squares) and L^p/a = 50 (black circles). 
The red continuous lines correspond to the continuum limit, 



eqs. (761 and (77 1 



variable conjugated to the harmonic number n, there- 
fore the structure Ad-{6) =fct{a/V9L^) corresponds to 
A(T„ =ict{y/na/L^p). We now extract this function. Us- 
ing the path integral formulation it is straightforward to 
get the structure 



Aa{0) - Act(O) = -2 / dtV{x,x; t) 
Jo 



(80) 



^ieAf[Ct 



Ct 



1 e 



-t/Ll 



where A/'[Cf] is the winding number of the closed tra- 
jectory. At large times the return probability coincides 
with the one of a plane, V{x, x; t) ~ (appendix F 3 1 ; 



A 

sinh Xt ) 



-It _ 



it can be obtained from eq. (|76| thanks to an inverse 
Laplace transform. UsingS3 

^{\ + ^. 
corresponds to 



,ie7V[C; 



we deduce that expression (77 1 

9t 



1Ct 



A Fourier 



2a2 sinh(6lt/2a2 

transforirPS gives the distribution of the winding num 
ber, plotted on figure [17] 



Tra 



2t cosh^(7raWA) 



(81) 



We have recovered the well-known Levy law for the dis- 
tribution of the algebraic area A = No? enclosed by 
a planar Brownian motioEp^ZHZlZS Yov i » and 
n 3> 1, the return probability conditioned to wind n 
fluxes is therefore expected to behave as P„(x,a;;<) ~ 
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V{x,x;t)Qt{n) : 



(82) 



8^2 cosh^inna^/t) 
A Laplace transform gives the corresponding harmonics 



a 
■inn 



where 



•^3(0 



Q-i/v 
2 — 

cosh y 



We extract the following limiting behaviours 



(83) 

(84) 

(85) 
(86) 



The constant is estimated numerically : — 1-239. 
The tail of the distribution corresponds to 



Aa-r, 



(2n)3/4 Y ttL, 



/Srrn aj 



(87) 



for <^ ^/na. The saturation of the harmonics for 
— > oo is given by : 




for ^/na <C Li^. 

Harmonics reach a finite limit for — > cx). It 

is interesting to compare the MC of the chain and the 
MC of the square network. For the chain, the behaviour 
of the MC near zero flux Ad^'=^^'''\9) - l/^ + cste, 
eq. ( 55 1 , is related to a weak logarithmic divergence of 



the harmonics Aan^'^^'^^ ~ \n{L^/nL), eq. (61 1. For pla- 
nar networks the WL correction presents a weaker di- 
vergence at zero magnetic field Aa{6) ^ In 0. There- 
fore /g d6'ACT(6') < oo and the harmonics reach a fi- 
nite value in the limit oo. Let us compute this 
value. The singular behaviour Ad-{6) ~ Cste -I- ^ln|6'| 
near zero flux is expected to dominate the harmonics 
behaviour Act„ ~ J^^dO In 9 cos {n9). The typical scale 
over which In 9 varies is 9\nl/9. If we deflne 9c by 
9clnl/9c — i/n, we see that the integral is dominated 
by the interval [0,6*^]. We have 9(, ~ l/(nlnri), whence 

l/n. More precisely, we 
A5-„ ~ -j^ for n> 1. 



"1/(71 In n) 
p 

have obtained above 



Numerical calculations.— The MC is computed nu- 
merically as a function of the reduced flux 9 for rational 
fluxes 9 — 2'Kp/q. As recalled in appendix G 3 the compu- 
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Figure 13: Numerical calculation of ten first MC harmonics 
A(T„ for a square network, computed for L^/a between 0.5 
and 50. 



of a tight binding Hamiltonian on a square lattice sub- 
mitted to a magnetic field, the so-called Hofstadter prob- 
lem. For rational flux 9 = 2Trp/q, this spectrum presents 
q bands determined by the polynomials Pp^g(e). For ex- 
ample, band edges correspond to roots of Pp^q{e) = ±4. 
For small q (in practice we choose ^ 8) the MC is com- 
puted by using eq. ( 75 1 . For large q (large number of 



Hofstadter bands)^, we use a more efficient procedure and 
rather follow Ref.^^ : we neglect the dispersion of Hof- 
stadter bands, according to which eq. (GIO) reduces to 

TT^Ny I Nijfip,,) } - I Sr=l 4cosh(y7a)-|-£, ' ^here Er 

designates the position of the band. The weak localiza- 
tion correction is represented on figure [12] as a function 
of the reduced flux 9 for three values of the ratio L^p/a. 
As this latter increases, the MC becomes sharper around 
zero flux, according to the above discussion and the har- 
monic content becomes richer. The MC is computed in 
this way for different values of the phase coherence length 
ranging from L^/a = 0.5 to 50. For each curve the first 
ten harmonics are extracted and plotted as a function of 



L^/a on figure 13 



In order to analyze the numerical results, we use 
the discussion of the above paragraph on the contin- 
uum limit. In the li mit L;^ ^ a we expect the scaling 
rit ~ t/a^. On figure 14 we plot ^^A(t„ as a function 
of l/\/f = L^/ay/TTn (note that the scaling nt ~ t/a^ 
is only expected for > a when we reach a "two- 
dimensional limit" ; for ^ a we rather expect the 
scaling corresponding to the isolated ring nt ~ y/t/a). 

After rescaling, all curves of figu re [T3] collapse onto 
each other as we can see on figure |14| (at least in the 
domain > a). Some significant deviation from ex- 
pression ( 83 ) occurs only for n = I. In order to analyze 
the behaviour for largest more precisely, harmonics 
are re-plotted as funtions of the variable ^ln(fo/0 on 
the inset of figure 14 : we check the linear behaviour 



tation of the MC is related to the study of the spectrum 



with this variable. Surprisingly, the continuum limit can 
be considered as a very good approximation already for 
n ^ 2. 



17 




0.1 0.2 0.3 0.4 0.5 

I ine„/i) 



Figure 14: The harmonics represented on Fig. \13[ are plotted 
here as a function of the variable — L^/a^im. The ma- 
genta dashed line corresponds to the continuum limit : func- 
tion J-^i^,). Inset: Same functions as a function of (^ln{(,o/0 
with ^0 =i 1-239. 



Remark : Brownian path/random walk. - We have shown 
that the distribution of the number n G Z of cells en- 
closed by a Brownian path in the square lattice is very 
close from the Levy law describing the distribution of 
the algebraic area enclosed by a planar Brownian motion 
(continuum limit) already for n ^ 2. It is interesting to 
point out that this remark also holds for the number of 
cells enclosed by a discrete random walk jumping between 
the different nodes of the square 

latticePEf. 



B. Model B 

1. The two-dimensional limit 

The thin film.— Let us first recall so me kno wn results 
for the plane (or thin film of thickness 5^ 3 | 33 | 35 | ^..^^ 
mensions the diffuson presents a logarithmic behaviour. 
The function W behaves in the same way, with a cut - 
off at small scales at the thermal length Lr^MMMB ■ 
W{f,f') = ^ln(||7^-?^'||/LT) for \\f-r'\\ > Lt. There- 
fore the functional governing decoherence behaves as 



2e^T 1 , 
1-^7^^'^ 



Dt 



(89) 



We recognize the sheet resistance^" _/?□ = 1/(6(To) of the 

film of thickness h. The phase coherence (Nyquist) time 

is evaluate d from V\Ct \ t ^ \. We obtain the temperature 
dependenc^^l^lM] 



1 

—film 



T In 



2Ru 



(90) 



valid for Lj- ^ h. This behaviour was observed exper- 
imentally for thin metallic film'*^ and two-dimensional 



electron gaj^. We recall that 2d magnetoconductance is 
given bjffllllMl : 



{B) = 



2e^^ 
IT 27r6 



-In 



+ cste (91) 



where rg = (j>o/{8TTDB) and ^p{z) is the Digamma 
functiorPS (the additional factor 1/2 in the Digamma 
function,^ compared to ( [TT] ), is explained in ap- 
pendix G 2 1 . We may simply write Act ~ 
1 



h 2-Kb 



ln[min (r^''", Tg) /re]-|-cste. The small time cut- 
off Tf, in eq. (|9l| is introduced by hand to account for 
the fact that the diffusion approximation only holds for 
timesPlt > Tp. 



The square network.— For large time scale {T^p 



t > 



a ) and small magnetic fields (such that (f) <ti (po) the 
result for the network should coincide with the one for a 
plane. In this case the function W entering the decoher- 
ence rate is W{x, x') ~ ^ ln(||a: — x'||/a) where — 
is the distance between the two points of the network. 
The logarithmic behaviour is now cut off naturally at 
the scale a. Because the function W presents a smooth 
logarithmic behaviour, we extract the relevant time scale 
(phase coherence time) by following the same lines as for 
the plane. We write 




2e2r a , ^Dt 

= t — In (92) 

CTps 27r a 

where s = wb is the section of the wires of width w and 
thickness b (figure [T|) . From this expression we extract a 
time scale reminiscent of eq. ( 90 ) for the film : 



•N 



Dnct 

Rg 



'"Ly h/e 
2R? 



T In 



2 1 



(93) 



This result is valid for Lt ^ a. In the opposite limit 
Lt ^, a, the cutoff in the function W should rather b^^ 



Lt therefore 1/t^°* 



h/e^ 



T In 



r h/e 



However this 



latter regime seems less relevant from the experimental 
point of vie"wE^. The sheet resistance of the network is 



rjnct 



wbao 



SCTo 



(94) 



This characteristic time is reduced by a factor w/a, com- 
pared to the Nyquist time (90 1 obtained for a film of 
same thickness : ~ ^^n"^- We can also introduce a 
Nyquist length for the network L^* = ^/ Dt^P^ . related 
to the Nyquist length of the wire Ln hy 



1 

T net 




(95) 
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We expect that the MC presents the logarithmic be- 
haviour Act ~ ^ [In + Csn] which is cut off at very low 
magnetic field : Act ~ — Mn[min (L^*, Lg) /a], where 



Lb = V</'o/(47r6) is the 2d cutoff. 



2. MC Harmonics 

In the network, the diffuson behaves logarithmically at 
large distances Pd{x,x') ~ — ^ ln(||a; — a;'||/a). There- 
fore we expect that the relaxation of phase coherence is 
controlled by 



Vfit 



ttL% 



t ln{Vi/a) —f - — -g- t ln{LN/a) . 

(96) 

As for the plane we use the fact that the functional de- 
scribing decoherence weakly depends on trajectories since 
W{x, x') ~ In I |a; — x'l |. This suggests that the result for 
model B is given by performing, in the result for model 
A, the substitution 



1 



3a 
2^ 



T net 



(97) 



where the Nyquist length for the network is given by ( 95 1 
Using ( [87| , we get for the tail : 

'1/2 p-yg^a/L^V"' 



Act„ 



(98) 



(rlnl/r)^/'* g-ni/2a^/=^(Tl„l/T)i/2 



A similar substitution in eq. (|88|) leads to 
Act„ 



a 
Ann 



1 



\n\LN/a) 



(99) 



(100) 



to describe the saturation of the harmonics at large Lm /a 
(small temperature). 

We insist that since dO Aa{9) < oo the harmonics 
reach a limit for L^^ —^ oo (or Ljsf oo) which is indepen- 
dent on the decoherence mechanism. In other terms the 
magnetoconductance curve Aa{6) reaches a limit apart 
in a very narrow region of width 5 
zero flux (figure 12 1 . 



{a/L^Y around 




Figure 15: A metallic film is deposited on an insulating 
wire. This allows to study quantum transport in a long hol- 
low cyltndeS^. Two diffusive trajectories are represented with 
winding n — Q and n — I. 



Model A 



Let us first recall the well-known result for the weak 
localization correction computed within model jP^^. We 
denote by y e M the coordinate along the axis of the 
cylinder and x S [0, L] the coordinate in the perpendic- 
ular direction. The WL correction is written as a path 
integral over Brownian paths r^r) = (a;(r), 2/(t)) in the 
cylinder, where x{t) describes a Brownian path on the 
circle and ^(t) on M (figure 15 1 : 



Act 



-AAS 



-2 / dte"*/^^ 
Jo 

X r%r-(T)(5„,^[,(,)]e-^*'i^J^"' (101) 

Jr{a)=r 



dte"*/^*- 



-(nLy/4t 

Ant 



(102) 



VIII. THE HOLLOW CYLINDER 



where ^^0(2;) is a modified Bessel function. Therefore : 



We have noticed that a network made of a large num- 
ber of rings realizes disorder averaging. Another natural 
way to realize this averaging is to consider a long hollow 
cylinder of perimeter L (lon ger than L ^) submitted to 
a magnetic field along its axis^^IiSIiZllIl^ study below 
how the original result of AA^^ obtained within model A 
is modified when decoherence is dominated by electron- 
electron interaction. In this section, it is natural to define 
the reduced dimensionless conductivity ct as ct 
where b is the thickness of the metallic film. 



Act 



-AAS 



27TnL 



g-^ni/L^ for < \nL\ (103) 



In (L^/nL) 



for InLl < 



(104) 



These results are very s imilar to the one obtained for the 
chain of rings ( 60|61|62 ) . This is due to the similar wind- 
ing properties, what was already noticed after eq. (64 1. 
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Model B : e-e interaction 



We have now to consider the path integral 

(105) 



A5-„ = -2 / dt 
Jo 

X / 'Dr{T) 5n,J^[x{T)]Q 
Jr{0)=0 



where we have used translation invariance along the two 
perpendicular directions in order to deal with a path in- 
tegral with action local in time, in a similar way as for 
the ring : cq. ( Al I. 



1. The function W 

The cylinder is translation invariant in the two direc- 
tions, therefore we may write VF(r, r') = W{r — ^',0) 
with W{f,0) = Pd{fc,0) — Pd(^, 0), where fc is a short 
distance cutoff of order Lt (we will see that the direction 
of the vector plays no rolej^. 

In order to avoid the divergent contribution of the zero 
mode of the Laplace operator, we start by considering the 
solution oi {-f — A) P — S : 



1 f^°° 

^ rjj J — OO 



g2i7ma:/L+ifciy 

2^ 7H-(^)2 + fc2 



(106) 



2L 



E 



32i7rna;/L-^7+(?f=)2 \y\ 



Next we take the limit 7 ^ in 

W{r, 0) = lim \P{r^, 0) - P(f, 0)] 

7^0 



2L ' 2tt 



1 °° 1 

-T- 



2'!mXc _2lL2jJc 

COS — - — e ^ 
n \ L 

- cos — - — e ^ 

Lj 



(107) 



(108) 



(109) 



We finally obtain 



In 



1 -e-T 



2tiLt/L 



(110) 



where we used that \\rc\ \ = Lt <C L. We can check that 
this expression reproduces known results in two limits : 
for \y\ ^ L, we recover the Id form W ~ \y\/{2L). For 
|y| < L we obtain the 2d result W ~ ^ ln(||^|/LT). 



2. Harmonics 



The first term of eq. (110) originates from the Id mo- 



tion along the cylinder. To this Id motion, we can asso- 
ciate a Id Nyquist time similar to the one obtained for 
the wire, eq. (Isl, 




(111) 



that coincides with ([5| in which the section is taken 
as s = bL. 

We consider first the high temperature limit L 3> 
LnacL = ^jDTN.id = (^^^^)^/^. We remark that for the 
harmonic rt = 0, trajectories very unlikely wind around 
the cyhndcr and we can use W ~ ^ ln(| |r| |/Lt) for 
y <^ L. Therefore the calculation of the path integral 
( 105 1 corresponds to the one for the film, Actq — Acr*^'™, 



with the time scale ( 90 1 



Next we consider non zero harmonics n 7^ 0. In this 
case the trajectories have a small extension along the 
"^ire \y\ < L^^id <^ L and we can neglect the \y\ in the 
exponential in eq. ( |110| ) (see figure 16 1. Therefore we 
perform the substitution : 



with 



(112) 



(113) 



This approximation allows us to factorizc the path inte- 
gral as : 



A5-„ ~ -2 



dt 



rx{t)=0 



'x{0)=0 



rvW=o 



'y{o)=o 



Vy{T) e 



(114) 
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The first path integral runs over trajectories encirchng 
the cyhnder. Therefore we can replace W{x) by its aver- 
age /g ^W{x). This approximation is justified by the 

fact that W has only a logarithmic dependence. This 
simplify the calculation by substituting the functional by 
a constant : 



'o 



t 



cyl 



(115) 



where we have introduced the time scale 



cyl 



h/e- 



■Tin 



(116) 



This time is reminiscent of the Nyquist time for the film, 
eq. ( 90 1 , but the two times differ by the argument of the 



logarithm. 

The second path integral precisely coincides with the 



one for a wire 



Act„ ~ -2 



/47rt 



'.Ct given by (24 1. Finally 




1" m|t/Tw,ld 



(117) 



which leads to the series (for n^Q) 




(118) 



where the times are defined as 

1 dcf , 1 



tnm 



cyl 



(119) 



(we recall that Wm's are zeros of Airy function Ai'). 
This expression assumes that L 3> L^^id- We show 



that (1181 it also valid for the other regime L <^ LN,id 



in this case the path integral runs over trajectories such 
that \y\ ^ L (see figure 16 1, therefore, in (110 1, the ex- 
ponential damping suppresses the x and y dependence in 
the logarithmic of W what leads to the same conclusion 
for the two regimes since W + ln(L/27rLT)- 

In order to analyze the two limiting cases into more 
details it is convenient to relate the two times as : 



Tnm 

cyl 



L 



L 



N,ld 
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L 



(120) 



Since the two lengths i^2d are related, we 

just have to consider two different regimes. As it is clear 
from eq. (119), the harmonics are always controlled by 



the smallest scale among 2d ^'^^ ^ 



• High temperature L^^id <C L (then i^2d L^^id)- 
The WL correction is dominated by non winding trajec- 
tories in this case 



27r 



film 
N 



cste 



(121) 



involves the 2d Nyquist time ( 90 1 . Considering the os 



dilating part of the MC, only the first term of the series 
dominates. The harmonics are governed by the smallest 
length among i^2£i ^w.id • 



Ao-„ 



" cyl 
^N.2d 



2\ui\L 



N,ld 



I L (T In T 



,1/2 



(122) 

for LM.id ^ L. Note that this result is reminiscent of the 
result (27 1 for a ring Aa™^ ~ q-uL^/'^t^/^ . gome 
logarithmic correction it presents a similar T^/^ in the 
exponential for the similar reason (related to potential 
fluctuations seen by winding trajectories). However the 
L dependence differs from the one of the ring. 

• Low temperature L ^ LN,id (then Ljv.irf ^ ^'N^2d)~ 
The harmonics involve Ljv.id, the smallest length among 



L^2rf ^'Hd Liquid- Eq. (118) coincides with the one ob- 



tained for the chain of rings (68) : 

Act„ ~ -2F2(7iL/Ljvad) 



(123) 



For intermediate coherence length, L <^ i^v.id ^ nL, 
eq. (118) gives 



Aa„ 



1 



|ui|3/2 



-|tii|i/^riL/Ljv,ni 



-nLT 



1/3 



(124) 



For largest phase coherence length nL <C iw^id, we 
may use the form derived in section |VI| 



(125) 
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^NAd 

nL 



cyl 



where the constant, Ccyi 
tionlVll 



0.51, was introduced in sec- 



Discussion.— It is worth emphasizing the similarity be- 
tween the results for the cylinder and for the networks. 

In section IIVI we have seen that the MC harmonics of 
a weakly coherent ring probe two length scales L^ (x 
y-i/3 or Lc oc T~^/^. For the lowest temperatures the 
surrounding network matters and another length scale 
emerges : the MC of the square network involves a unique 
time scale 1 /t^°^ ~ Tin 1 /T, eq. ( 95 1 , reminiscent of the 



2d Nyquist time ([901. 

For a cylinder the MC also probes several time 
scales. At high temperature the zero harmonic related 
to non-winding trajectories probes the 2d Nyquist time 
~ ^Je^ ln[/i/e^2i?n] whereas the nonzero harmon- 

^Tlii[L^/Lj.]. The main 



ics probe the time 



cyl 
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Figure 16: Trajectories contributing to first harmonic for 
"high temperature" (left) and "low temperature" (right). 



dependence of the corresponding length L'^^^^d ^ T~^/^ 
has the same origin as for a single ring and reflects that 
winding trajectories feel fluctuations of the potential over 
length scale given by the perimeter (figure [T6| left). For 
lower temperature, trajectories diffuse along the cyhnder 
over length scale much larger than the perimeter and the 
WL correction is controlled by a unique length L^^id, cor- 
responding to the usual Id Nyquist time rjv,i(i oc T~^^^. 



IX. CONCLUSION 

We have studied the weak localization correction in 
metallic networks and in a hollow cylinder. This study 
relies on a detailed analysis of the winding properties of 
closed Brownian trajectories in these systems. We now 
summarize our results. 

We first recall the behaviour of the probability to re- 
turn to the starting point after a time t for trajecto- 
ries conditioned to wind n rings. In the short time limit 
t < L^, we have 'P„{x,x;t) ~ p„-^e"("-^)^/^* where 
Pn depends on the network : pn ~ I for the isolated ring, 
Pn — (I)" " for the ring connected to Na long wires 

and Pn — ^2"+"iwf chain of rings. For the square 

network, there is no close expression but a systematic ex- 
pansion may be found in Ref.^i^ with the trace formula 
of Ref.EH. 

Network V{x,x;t) Qt{n) q(x) 



(47rt)l 



xxxx 



1 



nL \ 
s/2i) 



- 1 1-2 











M- a 2\/Zo?- /'n2\/3a2^ tt 1 










4irt t ^\ t I 4\/3 coBh2(^^) 



Table I: Return probability and distribution Qt{n) — ^^^'^.'^^ 
of the winding number in the large time limit t ^ , The 
function is defined m eq. (41 1 



At large times t ^ the typical winding number 
scales as rit ^ (t/L^)", where a is a network dependent 



exponent. Introducing the return probability V{x, x; t) = 
X^n "Pnix, x; t), we may write the winding probability as 



T^n{x^ x, 



V{x, x; t) 

WWW 



n 

C2 (t/L^)^ 



(126) 



where J dxq(^x) — 1. The dimensionless number C2 en- 
sures that / dxx^q{x) — 1. Since V{x,x;t) ~ t^^^^ 
where d is the effective dimensionality of the network, 
we may also write 



"^n {x^ X, t) 



1 



ta+d/2 ■='Vt' 



(127) 



(L may be re-introduced by dimensional analysis). The 
function q(x) is given for the various networks in the 
table |l] and represented on figure 17 Surprinsingly the 
functions for the connected ring and the plane are very 
close ; they only differ in the wings when functions are 
exponentially small. 




Figure 17: Function q{x) for the connected ring (black con- 
tinuous line), the chain of rings (red dotted line) and square 
network (blue dashed line) in semilog scale. Inset : Same 
functions in linear scale. 



We have analyzed in details the harmonics of the mag- 
netoconductance oscillations obtained when decoherence 
is described by a simple exponential relaxation {model 
A). In the limit of large coherence length compared to 
the perimeter of the rings, the scaling of the harmon- 
ics can be easily understood from the Laplace transform 
Act„ ~ dtVn{x,x;t) e~*^^f . We see that the time 
scale coincides with t ^ L'^. We deduce from ( 127 1 that 
harmonics are of the form 



1 



(128) 



where <i>(a;) is a dimensionless network dependent func- 
tion (the perimeter L is easily reintroduced by reminding 
that A(T has dimension of a length) . We may then sum- 
marize for each geometry : 
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• For the isolated ring {d = 0, a = 1/2), the form of 
the harmonics Adn L^^i ,-(nL/ L^) is related to 
the scaling nt ~ t^^^. 

• For the connected ring {d = 1, a = 1/4) : 
AfT„ ~ ^ L^L <^c.v{nyj L/ Lip) can be understood 
from rit ^ t^^^. 

• For the chain of rings (d = 1, a = 1/2) : Aa-„ ^ 
^chain{nL/Lp) reflects rit ^ t^/"^. 

• For the square network (d = 2, a = 1) A(t„ ~ 



^ <&s.n(f^i^/i^) originates from rit ^ t (here the 
harmonics were not written exactly under the form 
(128), but in terms of the function <b{x) — x^{x) 



in order to emphasize that harmonics reach a finite 



value for 



oo). 



The precise behaviours for the harmonics Ag„ are sum- 
marized in tableau (we recall that A(7„ ~ AtT„ apart for 
the connected ring for which Agn ~ y^A(T„ where la is 
the length of the connecting wires). 



Model A (exp. relax.) 



Model B (e-e inter. 



Regime <C L 



an 



^l/2g-nL/L, 



for L"^' < Ljv,id < L 



Regime ^ L 



r3/2 



for < L^/L L^^^ 



for < Liv /L 



for > Ly,/L L^/^e~"""V^/^" for > Ljy/L 



foTn<^Lip/L ln(I/jv/nL) 
for n > L^/L e-"'^"'^/-^" 



for n <^ Ln / L 
for n ^ Ln / L 



for ^ L^/L 
for 01 S> L^/L 



idem for -+ L^' 



idem for L^°' 



or 



ln(Z/^/nL) 



for n <: I/,p/L ln(Ljvad/'^^) 
for n > L^/L e"'"""'^''-^"." 



for n < Ljv,id/i 
for n > Ljv,id/L 



Table II: Harmonics of MC Ag„/g for different networks. 
In the high temperature regime <^ L, winding trajec- 
tories cannot explore more than a single ring and harmon- 
ics do not depend on the network. Dimensionless constants 
are m = ~ 1.234, K2 = V2\ui\^/^ ~ 1.421, K.3 = 

2"^''^|tti|^/^ ~ 0.801 and K4 = M^^'^ ^ 1-009. The vari- 
ous Nyquist lengths are Ln ~ T^^^^, LJ,"' ~ (Tlnl/T)"^''^ 
anrf ~ (rinr)-^''^ 



For each situation we have also discussed the effect of 
decoherence due to electron-electron interaction {model 
B) , the dominant phase breaking mechanism at low tem- 
perature. As recalled at the begining of the paper, this 
mechanism requires a refined description : the simple 
exponential decay of phase coherence is replaced by a 
functional of the trajectories, eqs. ( 15|18 ). In networks 
of quasi- Id wires the decoherence due to e-e interaction 
is controlled by the Nyquist length Ljv cx; T~^/'^ . 

In the "high temperature" limit Lm <C L the fact 
that trajectories with finite winding number and tra- 
jectories with winding n = do probe different length 



scales is responsible for the emergence of two length 
scales Ln cx T^^/^ and L^ oc T^i/^ (q^. j^shn ^ y-1/2 



and L 



cyl 
N.2d 



(Tlnr)-i/2 for the cylinder). The mod- 
els A & B give different dependences in the phase co- 
herence length : Atrf^' -L^e~'^^^^f and Acr^^' 

— Ljve^"'-'^/^"-'^'' The exponential decay of harmonics is 
almost independent on the network. 

In the "low temperature" limit Ljv 3> L, all trajec- 
tories probe the same typical scale, irrespectively of the 
winding. However this length scale depends on the ge- 
ometry : Ltv oc T~^^^ for the chains of rings and the 
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hollow cylinder, and L^^* = {^^Y'^\n~^''^{LN /a) - 
(Tlnl/T)-i/2 foj. ^-^^ 

square network. As a function of 
the phase coherence length, models A & B predict har- 
monics of similar form strongly network dependent. We 
have compared harmonics as a function of the phase co- 



herence length for the different networks on figure 18 (for 
model A). 

All results are summarized in table|llj We have plotted 
the WL correction to conductances for the three different 
networks on figure [TSj 

An experimental verification of these predictions would 
be interesting and would confirm our understanding of 
decoherence due to electron-electron interaction in com- 
plex geometries. In particular an interesting and clear 
experimental test would be to compare the MC harmon- 
ics for the chain of rings for independent rings and coher- 
ent rings (networks of figure |9| in the "low temperature" 
limit Ln ^ L. The experimental analysis of the MC har- 
monics for the square network in this limit seems more 
difficult due to the fact that harmonics reach a value in- 
dependent on the phase breaking mechanism. Therefore, 
contrarily to the chains of rings, the MC harmonics for 
the square network are less sensitive to the model of de- 
coherence for large phase coherence length. 



I 































































































xlU 






















LJL 



00 

I 

































































































xlO 


•< 

















Figure 18: Comparison between harmonics Agn/g for differ- 
ent networks in the low temperature regime > L with n = 1 
and n — 5. For the connected ring (continuous black line) we 
choose la/L = 10. The blue dashed line corresponds to the 
chain of rings. The result for the square network (dotted red 
line) has been multiplied by a factor 10 for visibility. 



Acknowledgements 

We thank Christopher Bauer le, Helene Bouchiat, 
Markus Biittiker, Richard Deblock, Jean Desbois, Meydi 
Ferrier, Sophie Gueron, Alberto Rosso and Laurent Sam- 
inadayar for stimulating discussions. 



Appendix A: A USEFUL PROPERTY OF 
WINDING BROWNIAN TRAJECTORIES 



The difficulty for computing the path integral ( 19 1 lies 
in the time nonlocality of the action. In this appendix we 
show how it is possible to get rid of t ime nonlocality in 
certain cases, as explained in Refs.l^lSli. For that purpose 
we demonstrate the identity 



x{t)=0 
x(0)=0 



Mb 



e 2 Jo 



(Al) 



where x(t) is a Brownian path on the circle (here identi- 
fied with the interval [0, 1]). The identity is valid for any 
symmetric and periodic function : V{—x) = V{x) and 
V{x + n) = V{x) for n G Z. 

Demonstration for n = was given in Ref.'^'^I!^ where 
we pointed that, for a Brownian bridge on M (a;(r), ^ 
r ^ t\x{0) = x{t) — 0), we have the following equality 
in la^^ : 

x{t) ~ x{t - t) ^ = \{2t) (A2) 

The proof lies on the fact that we can relate the bridge 
to a free Brownian motion (Wiener process) {W{t), ^ 

T < t, W{0) = 0) : x(t) * = ' W{t) - ^W(t). 

Here we generalize this relation when x{t) lives on the 
circle [0, 1] and when we constraint the winding number. 
Let us unfold the ring in order to work on R. A close 
path winding n times around the ring is related to the 
following path living on the real axis : (a;„(T), ^ r ^ 
t\ Xn{0) = ; Xn{t) = n) that can be writen as 



Xn{r) ''='w{T) + -^{n-W{t)) ' = 



Xo{T) + n- (A3) 



Xq{t) is the Brownian bridge. It is now easy to show 
thaiPl 

Xn{T)~Xn{t-T) Xn{2T)-n for T < t/2. (A4) 

Since x„(r) — Xn{t — t) is argument of the periodic func- 
tion, the integer shift can be forgotten. The symmetry 
V{x) = V{—x) ensures the equality of contributions of 

intervals jj^^ and J^^^. It follows that 



/ dTV{x^{T)-X^{t-T))"=' f dTF(x„(T)) 

Jo Jo 



(A5) 
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which demonstrates eq. (All 



Infinite wire : Using (Al I, we see that the path integral 



( 20 1 involves an action local in time 



1 r°° 
^ Jo 



X / 'Dx{T)e 



(A6) 



that can now be computed. We obtain Pc{0, 0) = 
2 a'/'I-yL^^ ) derived in RefP (numerical factors are 
incorrect in this reference). 



Isolated ring : The function W{x, x') is given by eq. (E9 1 
The path integral (19) can be rewritten as 



die-''* 



x{t)=Q 

Vx (5^ 

x{0)=0 



-/od^[^+i#(i-¥)] 



(A7) 



that can be expressed in terms of Hermite function^. 




Figure 19: Appropriate contour deformation in order to esti- 
mate (Bll. The dashed area corresponds to the region where 
Re[(p{zJ]< 0. 



(note that a factor 1/2 is missing in 

Finally the relation to the function q{x) introduced in 
the conclusion requires the two integrals Jp°°d^\E'(^) = 
l/^/^and/„°°d^^2^(0-2. 



Appendix B: THE FUNCTION *(^) 



We analyze several properties of the function ( 41 1 , that 
we rewrite 



e 'T 



(Bl) 



where Lp{z) = z'^ + 426"'''/'' and A = (C/4)'^/^. The value 
of the function at the origin is ^(0) = ~/|^- 

The asymptotic behaviour for ^ ^ 1 may be stud- 
ied by the steepest descent method. <p'{z) — has 
three solutions z„ = e''"/4+2i»'^/3^ ^ith n = 0, 1, 2 (fig- 
ure [l9|. An appropriate contour deformation in the com- 
plex plane of the variable z must remain in the region 
where Re[(p(z)] > 0. This domain can be easily de- 
termined by performing a rotation z — wG^^^'^ : writ- 
ing w = u + iv we have Re[(y5(z)] = — + 6u^v^ — 

V + Au that vanishes for v = ± 
The domain where Re[iy9(z)] > is represented on fig- 
ure |19[ This shows that the contour can only visit 
zq — e'^/^. The integration over M+ is replaced by 
integration over the segment A from the origin to zq 
and the cont our C issuing from zq and going to in- 
finity (figure 19). Noticing that dzz^ e''^'^^^) is 



purely imaginary, we are left with the contribution of 
the contour C only. We now use the steepest de- 



scent method ^{C) 
^A3/4Re[e-'Tl 



^ A3/4 Re[e-'T / dz z2 e-'^ '^(^)] 



_2|_g-A¥>(^o)]^ where the 1/2 is due 

to the fact that the contour issues from the stationary 
point, hence 



» 1) 



l/3p-3(?/4)''/=' 



(B2) 



Appendix C: HYPERGEOMETRIC FUNCTION 
F(|,n + |;n + l;|) 

This appendix is devoted to the study of the hyperge- 
ometric function _F(i,n-|-^;n-|-l;^). Our starting point 
is the integral representatioiJ^ 



,1 1 



At 



2y^U'^^ 



^{l~t){l-£,t) 



(CI) 



We recall that the hypergeometric function is regular at 
the origin F{a, /3; 7; 0) = 1. Note that the Euler /3 funtion 
_B(i,n+ i) — Y^r(n-f \)/n\ is well approximated by 

i3(|, n -f ^) ~ \J 1^ jn in the large n limit. 

In order to analyze the behaviour of the hypergeomet- 
ric function for ^ — > 1 we rewrite the integral of eq. (CI I 
as 



VK^-^ + ^t) 



(C2) 



1 



f 



dt 



1 



1 , (l-t)"-l/2_l 

d< / 

VW+et 



where e = 1/C ^ 1- The first integral is 2argsh(l/y^) — 
ln(4/e) + 0{e). The second integral reaches a finite limit 
for e — > 0, expressed in terms of the Digamma functioiJ^ 

Jo "" i^u~^ — "0(1) — V-'('T' + 5)- We can show that 
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correction to this constant is of order £ln(£), therefore 
^ • -t- 1 • f I = ^ 



4^1 i?(in+i) 



(C3) 



In 



The behaviour (C3) only holds for n not too large, n <^ 
(1 — 0^^- opposite case n » (1 — £_)~ ^ ^ 1, 

the factor (1 - t)"-i/2 ~ e-("-i/2)t selects 
an interval of width l/rt ^ (1 — ^) and we can neglect 
the quadratic term below the square root. Therefore 

__Q-in-i/2)t^ Finally 



eq. (|C2]) is ~ dt 



1 



(C4) 



for 77 > (l-^)-i > 1. 




Figure 20: Comparison between the two sides of equation (C6l 
for n = 1, 2, 10 and 100. 



We now prove a useful relation between the hyperge- 
ometric function and the MacDonald function (modified 
Bess el function). If ^ = ~ 1 — 2z/n, the integral 

( C2 1 may be rewritten as 



" dt (1 
n 



1-1/2 



n \ n 



dt , . (C5) 

^/t{2z + 1) 



We recognize an integral representation of the MacDon- 
ald function^''. Therefore, for z ^ n, we can write 



,1 1\ /I 1 

B\ -,n+ - F\ -,n+ -:n- 
2' 2 I 2' 2' 



l;e- 



(C6) 

The r.h.s. describes the crossover between (C3) and 



(C4|. We compare the two sides of equation (C6) for 



different values of n on figure [20| 



Appendix D: LAPLACE EQUATION IN 
NETWORKS : SPECTRAL DETERMINANT 



In this appendix we introduce an important tool, the 
spectral determinant, used to study some properties of 



the equation 



(7- A)P(x,a;') = 5{x-x') 



(Dl) 



in networks. 

The spectral determinant is formally defined as >S'(7) = 
det(7 — A) — 011(7 + '^") where {En} is the spectrum of 
the Laplace operator —A (in the presence of a magnetic 
field, A ^ [V — 2ieA{x)]'^). Despite this operator acts 
in a space of infinite dimension, the spectral determinant 
can be related to the determinant of a finite size matrix, 
of dimension equal to the number of vertices. This ma- 
trix encodes all informations on the network (topology, 
lengths of the wires, magnetic field, boundary conditions 
describing connections to reservoirs). Let us label ver- 
tices with greek letters. Iaf3 designates the length of the 
wire (a/3) and Oa/s the circulation of the vector potential 
along the wire. The topology is encoded in the adjacency 
matrix : aa/3 = 1 if a and /3 are linked by a wire, aap = 
otherwise. We consider the case where Laplace opera- 
tor acts on functions ip{x) (i) continuous at the vertices 
satisfying (ii) J^/s '^af^V^a/si^) = where ipapix) des- 

ignates the component of the function on the wire (a/?) 
and ipa its value at the vertex. Self-adjointness of the 
Laplace operator is en sured if Aq S M (more details may 
be found in Refs.'^^'Sil). Aq, = 00 corresponds to Dirich- 
let boundary condition at the vertex and describe the 
case where a touches a reservoir through which current 
is injected in the network. Aq = for internal vertices. 
The interest of mixed boundary conditions (finite Aq) is 
illustrated in appendix [F| We introduce the matrix 

Maf3 = Safi ^Aa -|- ^7 X! ^7^^ 



sinh ^laSi 



(D2) 



where the co nstrain s the sum to run over neighbour- 
ing vertices. TheiP^'^H 



n 

(a/3) 



sinh ~^la0 
V7 



detTW 



(D3) 



where the product runs over all wires. Despite the spec- 
tral determinant encodes the spectral information, it 
is also possible to extract some local information, like 
P{x,x), by small modifications of the matrix. This has 
been used in Ref.^^and is briefiy discussed in appendix [f| 
It is useful to remark that the matrix M can be used 
to express P{x,x') when x and x' coincides with nodes 
(this is always possible to introduce a vertex anywhere 
without changing the properties of the network) : 



P(a,/3) = {M- 



(D4) 



WL correction in regular networks. - In large regular net- 
works connected in such a way that currents are uni- 
formly distributed in the wires, we can assume that 



26 



weights attributed to the wires of the networks in eq. ([s]) /3 as TZa/B = {Va — Vp)/!. Therefore : 
are equal. In this case, a uniform integration of the 
Cooperon Pc{x,x) — {x\::p-^\x) in the network leads 
to a meaningful quantity (relevant experimentally) . The ctq^ 
Cooperon integrated uniformly is directly related to the 
spectral determinantP^I^ 



(A^o')«c« + (Mo')/3/3 



(E4) 



in 



dx /S.cf{x) 



network 



-2|-ln^(,) 



(D5) 



This equation provides a very efficient way for calculating 
the WL correction in arbitrary networks, when uniform 
integration of Cooperon is justified. 

WL correction in arbitrary networks. - In the most gen- 
eral case, eq. ([S]) requires to construct the Cooperon in 
each wire. A general expression was provided irP^ how- 
ever it is useful to notice that Pc{x,x) can also be ob- 
tained from a spectral determinant for a modif ied b ound- 
ary condition at point x. It was shown in Refs.l^Eilthat if 
we introduce mixed boundary conditions with a param- 
eter \x at X, then 



Pc{x,x) 



d 
OK 



lnS'(^-)(7)| 



(D6) 



Appendix E: CLASSICAL 
RESISTANCE/CONDUCTANCE 

We calculate the resistance between two vertices of an 
arbitrary network. We consider a network of wires of 
lengths with same sections s. In this case the con- 
ductance of the wire iiv is given by cros/lfiu- We introduce 
the matrix 



(El) 



whose matrix elements coincide with the conductances of 
the wires (up to the factor ctos). This matrix coincides 



ten and moreover with = 0, V /i for 



with the matrix (D2) if all fluxes are set to zero and the 
limit 7 — !■ is ta 
an isolated network. 

We now consider the situation where we inject a cur- 
rent at the vertex a. This current exits at vertex (3 (see 
figure 21 ). If we denote by the potential at /x, Kirch- 
hoff law at vertex /i takes the form 



CTqS 



Potential is therefore given by : 
/ 



(E2) 



(E3) 



Note that the matrix Aio is not inversible ; it is ex- 
plained below how to give a precise meaning to this ex- 
pression. We define the resistance between points a and 



Using eq. ( D4 1 we see that we can express the resistance 



in terms of the solution Pd of the equation — AP^ — S 



TZ{x,x') 



2 



Pd{x,x) + Pd{x',x') 



W{x,x') 



Pd{x,x') 



(E5) 



This demonstrates that the function W{x^x') defined by 
( |17| is indeed the equivalent resistance between points x 
and x' . 




Figure 21: Injection of current in a network (here a regular 
square network). 



Remark : M.o is not inversible.- It is easy to check that 



E(-^o)m. = 0, 



(E6) 



Kernel of the matrix is the vector (1,1,-- - ,1). Physically 
eq. (E6l ensures : (i) that sum of all currents arriving at 



vertex v is zero, (ii) currents are zero if all potentials 
are equal (equilibrium). This problem can be overcome 
easily by noticing that only differences of inverse matrix 
elements have appeared, in eqs. (E3l and (|E4|). We can 



always inverse the matrix M.q in the space orthogonal to 
the vector (1, 1, • • • ,1) and comput e difference s of such 
matrix elements. This is how eqs. ( E3|E4|E5 1 must be 
understood. 

In practice, an easier way to compute such differences 
is to regularize the calculation by computing the inverse 
of matrix for finite 7 (or at least one finite Aq.) and 
take the limit 7 — > (or Aq. 0) after having computed 
the difference of inverse matrix elements : 



7^0 



M 



(E7) 



Note that detM ^0 for 7 e R+* since Spec(A) c . 

This point is related to the fact that was already men- 
tioned in the continuum limit in order to construct the 
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function W{x, x') (in section VIII B 1 or in appendix [f|) : 
eq. (E7l is the analogue of eq. (108). In the continuum 



this problem is related to the fact that the Laplace oper- 
ator is not inversible in the space of functions satisfying 
Neumann boundary conditions corresponding to an iso- 
lated conductor. 

Example : Junction W in an isolated ring. - The relation 
between the function W and the resistance may be used 
in order to construct easily W . Let us consider the case 
of a ring of perimeter L. When the ring is connected at 
two wires at x and x' , the resistance TZ{x, x') corresponds 
to the one of two wires of lengths |x — and L—\x — x'\ 
put in parallel. We straightforwardly recover the function 
obtained in Ref.^ : 



Wnng{.X,x') 



1 



la;— L—\x~-x'\ 



1 - 



\x-x'\ 
L 



(E8) 
(E9) 



Note that if we consider a ring connected at reservoirs 
through arms of finite length, the function W{x^x') in- 
side the ring is modified^, since its construction needs to 
consider a four terminal device. In the limit of infinitely 
long arms we can neglect currents flowing through the 
leads and we recover the result of the isolated ring (ap- 
pendix of Ref.l^ . 



added to the matrix A^o, that now become inversible. 
These parameters describe Dirichlet conditions for the 
diffusion equation and permit to inverse the Laplace op- 
erator. The transport through the network is character- 
ized_by a conductance matrix whose elements are given 
b^ 



Qa'P' 



laa'l/3/3' 



(ElO) 



At flrst sight diffuson Pd{a,(3) {i.e. inverse matrix ele- 
ment (A^o ^)q/3) is related to both the conductance and 
the resistance. Note however that the matri x Aip (and 
the corresponding diffuson Pd{x,x')) in eq. (ElO I does 
account for the boundary conditions (Dirichl 
ary conditions at contacts, 
A^o (and Pd{x, x' 
network. 



Appendix F: SOLUTION OF THE DIFFUSION 
EQUATION IN SOME PARTICULAR 
NETWORKS 



et bound- 
primed vertices) whereas 
v')) in eqs. ( E4|E5 1 describe the isolated 



We consider the solution of the diffusion eq. (Dl I for 
the networks studied in this article. 



1. The ring with one or several arms 



IP' 



We consider a ring of perimeter L attached to an arm 
of length b connected to a reservoir (i.e. with Dirich- 
let boundary condition at its end) and submitted to a 
magnetic field. The spectral determinant i^ : 



Figure 22: A multiterminal network. Wavy lines represent 
contacts through which current is injected. Contacts corre- 
spond to vertices with primed labels. 

Classical conductance.- It is interesting to compare the 
formula obtained for the resistance with the one obtained 
for the conductance matrix of a multiterminal networlP^ 
We first s tress th at we consider two different situations. 
The eqs. ( E4|E5 ) give the potentials when the currents 
injected at vertex a and extracted at vertex /3 are fixed 
(figure 21 1. On the other hand, a conductance matrix 



allows to determine the currents through contacts when 
the potentials at these contacts are fixed (figure 22 1 . 

We consider a multiterminal network, connected to ex- 
ternal contacts through which current is injected at a set 
of vertices indicated with primes : a', /?',... The connec- 
tion is accounted for by introducing parameters Xa' = oo 



S{j) — —— sinh y^b [cosh y^Lcs — cos 6] 



(Fl) 



n 

where the effective perimeter is given by 
cosh y'^Lctf = cosh + - coth y'^b sinh y'^L . (F2) 

A systematic way for obtaining the spectral determinant 
of two subgraphs glued at one vertex from the spectral 
determinants of the subgraphs has been derived in " 
This allow to recover easily eq. (Fl|. 



• Introducing mixed boundary conditions at the 
node (vertex 0) we easily obtain 5'^'^"' (7) = 
^ sinh sinh ^6 + S'(7) from (|D3]). Using eq. ( [D6| , 



and performing a Fourier transform, we get the Cooperon 



pi") (0,0) = -iACT„(0) at the nodtP 



Pi"H0,0) = 



1 sinh y'^L 



2/7 sinh /7i, 



off 



(F3) 



In the weakly coherent limit we find Loff — L 



L^ln(3/2), whence Pi"^(0,0) 



1 

2^7 



(!) 
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Figure 23: A ring with one wire. We choose a Dirichlet bound- 
ary condition at the vertex 1 for simplicity (Xi =00). 



In the limit ^/^L <C 1 the effective perimeter is Lcff 



(F4) 



• In order to calculate the harmonics of the Cooperon 
in the arm (x is the distance from the ring, see fig- 
ure 23 a), we have to consider the spectral determinant 



for the graph with mixed boundary conditions at x, with 
parameter : 



5(^==)(7)-5(7) + 



sinh L 



sinh^(fe — a;) 
sinh b 



+ 2 sinh X sinh(& — x) (cosh Lcff ~ cos ( 



(F5) 



(for shorter notations we have omitted ^^7 in hyperbolic functions). From eq. (D6I we obtain : 



[x,x)-d„^o^ sinh ^6 sinh 776 J ^"'"^ 



(F6) 



for X e arm (figure 23 a). We recognize the first term as the result obtained for a wire of length b connected to 
reservoirs {i.e. with Dirichlet boundaries). In the limit b —^ 00 we have Pi"^(a;, x) ~ Pt\o, 0) e~2V7^ (for n^O). 

• If x is inside the ring (figure |23|b) the modified spectral determinant reads : 



S'(^-)(7) = 5(7) + — [sinha;sinh(L-x)cosh6 + sinh&sinhL] 

7 



and from eq. (D6) the Cooperon is therefore 



(F7) 



Pi") (a;, a;) 



1 + coth ^b 



sinh ^/7x sinh ^J^(Ij — x) 
sinh ■y/T'L 



Pi") (0,0) 



(F8) 



for X S ring (figure 23 b) 
If 



^ L we have Pc"'\x,x) 



2^ 



]Pi"^ (0,0) for x < L/2. In the bulk (for x and 
L-x > Ly) wehavePi"^(a;,a;) ~ |pi")(0,0) (figurejii]), 
where the factor | corresponds to the ratio of coordina- 
tion numbers at and at x. 

In the opposite limit L <^ <^ b the Cooperon is 
homogeneous inside the ring, as expected. 



From one to Na arms. - In the regime <C L, we have 
seen that the presence of one arm is responsible for fac- 
tor (|)" originating from the n crossings of the vertex. 
We immediatly deduce that the Cooperon in the ring is 



Pt\x,x) 



1 (2^ 

In order to study the regime ^ L, instead of con- 
sidering the network of figure [7] we discuss the case where 




Figure 24: x dependence of the Cooperon Pc"\x,x) in the 
connected ring for <^ L. 



all arms are attached at the same point in the ring. The 
calculation is more simple in this case. The two situ- 
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ations were studied in Reff^ where it was shown that 
the Cooperons for the two networks only sHghtly dif- 
fer in the regime <C L and are equal in the regime 
^ L of interest now. The effective length is now given 
by cosh ^J^LaS = cosh ^J^L + sinh ^J^L coth ^ b. 
The structure (|F3| still holds. When L <t: <t: b 

we find L^s — J~^^'^\/NaL, therefore Pc^\x,x) ~ 
\yj L^L/Na e~"V^°-^/-^'.', whose inverse Laplace trans- 



form leads to eq. (40 1. This Na dependence may be 



more simply obtained by noticing that, given the wind- 
ing probability Vn{x,x;t) for one arm, the one for Na 
arms is obtained thanks to the substitution n —>■ nNa 
and L L/Na- 



2. The necklace of rings 

We analyze the solution of the diffusion equation in 
a chain of Nr identical rings of perimeter L. Rings 
are attached in such a way that the two arms join- 
ing two vertices are symmetric. The chain is closed 
in order to form a necklace for simplicity ; as soon 



as the total length is smaller than L^, the results are 
insensitive to boundary conditions : periodic (isolated 
necklace) or Dirichlet (chain connected to external con- 
tacts). Let us label the vertices joining consecutive rings 
with greek letters a, P € {1, • • • ,Nr}. The matrix in- 
troduced in appendix |D] has the simple form Aiap = 
smh{^L/2) [Saf3 2cosh{y^^L/2)-aai3 cos(6'/2)] where 6* is 
the reduced flux per ring. With our convention the adja- 
cency matrix reads aaf3 — 5a.f3+i + Sa^f3-i. Its spectrum 
of eigenvalues is 2 cos{2mT / N,.) with n G {1, • • • ,Nr}, 
therefore the spectral determinant reads : 



Nr 



Y[ cosh(^i/2) - cos(6'/2) cos 



2Tm 



(F9) 



We replace the product by a sum by considering the log- 
arithm. The sum can then be computed in the limit 
Nr — !■ oo : 



^(7) = 

Nr^i 



2|cos(g/2)| sinh(V7L/2) ^argch( -t<^;,/;" ) 



Nr 



We now consider P{x, x') when arguments coincide with vertices : 



sinh ^L/2 



N,. 

E 



^r4:^/7 ^ cosh ^L/2 — COS 



27rn 



Nr 



-\a-f3\^L/2 



r 



(FIO) 



(Fll) 



OOGOC- 



-I- 



-I- 



-I- 



-I- 



X 



m L 

Figure 25: Infinite chain of rings. 

We introduce the notation x = {x, /) to locate points 
in the chain, where x is the coordinate along the axis (see 
figure) ; in the n-th ring x e [nL/2,{n+ l)L/2]. The 
discrete index / £ {u, d} indicates whether the point is 
in the upper branch or the lower branch. 

If X and x' are arbitrary positions we can consider sev- 
eral situations. 

• X and x' belong to different rings. The Cooperon is 



P{x,x') 



4^7 



(F12) 



that is half of the result for the infinite line. 



• X and x' belong to the same ring. 

P(i,i') ^ 



7^0 4^ 

We conclude that in time representation 



V{x,x';t) 



1 1 



'(x-x')^ /it 



(F13) 



(F14) 



(the relation is exact when x and x' belong to different 
rings). The 1/2 corresponds to the probability to end in 
one of the two arms of the rings. It ensures normalisa- 
tion : / dxV{x, x'; t) = d I ^xV{x, x'; t) = 1. 

The function W{x,x').— Expression (Fill can be 
used to construct Pd{x,x') by taking the limit 7—^0. 
We find for the function entering into description of de- 
phasing : 



1, 



W{x, x') — - \x — x'\ 



1 



(F15) 
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for X, x' in different rings, and 

W{x,x')=^-\x-x'\(l-^-^^]^W,,^^ (F16) 



for X, x' € same ring. In (F15I positions are measured 



along an axis along the chain ; the result is half of the re- 
sult for a wire. As it has been already noticed this result 



can be understood thanks to the relation (E5 1 : when the 
number of wires between each node is doubled, the resis- 



tance is divided by a factor of 2. In (F16 1 coordinates are 



relative to a unique axis inside the ring (it is understood 
that the expression is periodic) ; the result is exactly the 
result obtained for an isolated ring. Once again interpre- 
tation is easy : for an infinitely long necklace, when two 
external wires are plugged inside a ring, no current can 
flow out of the ring and resistance Ti{x, x') is not affected 
by the remaining rings. 



result call for two remarks. First, we see that the fact 
that the diffusion is constrained in the Id wires form- 
ing the square network leads to a continuum limit with a 
renormalized diffusion constant D* = 1/2, if the diffusion 
constant is Z? = 1 in the wire (more generally D* = 1/d 
where d is the effective dimensionality of the network). 
Second, the probability presents a prefactor 1/2 of simi- 
lar origin than the one in eq. (F14l, ensuring normaliza- 
tion condition. Each elementary plaquette, labelled by 
R = (n, m), contains 2 wires. Therefore each wire of the 
network can be labelled with R and an index / taking two 
values (for horizontal or vertical wires). We check that 
(F21l is correctly normalized : J^^^^^^^,^dxV{x, x' ;t) — 
EflE//wiro(i?,/)da;'^(a^,a;';i) ~ J2R'^aV{x,x';t) ~ 

2ay" - ii^e-^^(^--^')' ~ 1 

The generalization of ( |F18|F19|F20| ) to the other reg- 
ular planar networks can be done. We obtain : 



3. The square network 

We construct the solution of the diffusion equation in 
the square network. We use (D4| to express P{x,x') 
when the two coordinates coincide with nodes of the net- 
work. In this paragraph, nodes are labelled with a cou- 
ple of integers a — > (n, m) = R. The matrix A4 has the 
structure 



^/7 



sinh ^ 



/7a 



4cosh^adj^s,-aj^s, 



This matrix is easily inverted 
sinh y/^a 



R,R' 



(F17) 



(F18) 



2^/7 

dQ QiQ-iR-ii') 
2dB (27r)2 2 cosh y/^a — cos Qx — cos Qy 



where integral runs over wavevectors of the Brillouin 
zone. If we are interested in V{x, x'\ t) in the large time 
limit, t ^ a? , this corresponds to consider <C 1. In 
this regime the above integral is dominated by small Q 
and we get 

a f dQ eQ-iR-R') 



2 7o 27ri ^ ' 

where integral over Q has been extended to M^. We 
expect 'P{x,x';t) to vary smoothly on the scale of the 
lattice spacing, therefore we may write 



Vix, x'; t) 



a 1 
t^a^ 2 27ri 



where ||x — designates the distance in between 
the two points x and x' of the square network. This 



^ ' ' ^ t»a2 2 \ZJ 2TTt 



(F21) 



where the dimensionless factor ^ tan j , where z is the 
coordination number of the lattice, has the interpretation 
of the area of the Wigner-Seitz elementary cell Az{^ + 
^zfi), where Az is defined in the next appendix, divided 
by the number of bonds per cell. 



Appendix G: MC OF PLANAR REGULAR 
NETWORKS 



1. Planar regular networks 



We apply the formulae (D3) and (D5l to the case of 



infinite planar regular networks. Let us denote by z the 
coordination number of the network. The plane can be 
covered by only three different fillings : the triangular 
lattice (z = 6), the square lattice (z = 4) and the honey- 
comb lattice (z = 3). 

Let us mention few properties of regular fillings of flat 
surfaces : 

• the lattice of coordination number z is a tilling by 
regular polygones with p sides, p is related to the 
coordination number by^*^ (z — 2)(p — 2) = 4. 

• The area of a regular p-gone of side a is Az = 
a^\pcotg{-K/p) = a^^^^ tan(7r/z). 

• If boundary effects are neglected, the number of 
bonds B and the number of vertices V of the planar 
network are related by 2B = zV (z bonds issue 
from each vertex and each bond touch 2 vertices). 



(F20) The matrix (D2| takes the form M 



sinh a 



the matrix A* is given by 

Naf3 = SafS z cosh -^7 a — aafS e~'^°'' 



N where 



(Gl) 
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where reduced fluxes Ogfj de scribe a uniform magnetic 
field B. We use eqs. ( |D3|D5| ) to express the WL correc- 
tion A(7 = — ;^T^ ln S'(7), where the volume is related 
to the number of bonds Vol = Ba, 



2-1 

z 



( coth ^ 



1 



— L,„ smh — Tr < — 



(G2) 



6 = 47r</)/0o is the reduced flux per elementary plaquette. 
The computation of this expression requires the knowl- 
edge of the spectrum of the matrix N, i.e. of Oc^e"'^"'' 
(Hofstadter problem^"'). 



This result indeed coincides with the result for a plane, 



eq. ( 91 ), apart for a factor 2 in the argument of Digamma 



functions. We think that it is worth devoting a small 
paragraph to this point since it is related to some inter- 
esting property of the continuum limit of the result for 
the networks. 

The additional factor of 2 and the continuum limit.- As 
mentioned in section IVIII the WL correction in the net- 
work can be written with a path integral as 



/"OO 

Aa{e) = -2 / die-*/-^' 

"'0 



(G5) 



'0 

X Vx{t) ^' W[.(r)] 



2. Continuum limit 

We study the continuum limit : for ^ a the WL 
correction probes large scales. Moreover when the flux 
per cell is much smaller that the flux quantum, (p <ti (po, 
the result for the planar network coincide with the re- 
sult for a plane. Let us show that this is indeed the 
case. We note that the action of the adjacency matrix 
aafs on a smooth function can be replaced by the Lapla- 
1^ A + z. Therefore in the presence of 



'(V-2i eA)2 



cian : aa/s — > |/ 
a weak magnetic fleld Uafs e~'^°'^ 
z involves the covariant derivative. The spectrum of 
this operator is the Landau spectrum shifted by —z : 
Spec(-aa/3e-'^°") ~ {za'^eB{n+l/2)-z\n e N}, where 
each Landau level has a degeneracy cZll = ^ eKArea, 
"Area" being the total area occupied by the planar net- 
work. 

In the limit ^ a and 6^1, eq. (G2 1 is dominated 



by the trace which is itself dominated by the bottom of 
the spectrum of —Gap e~'^°'' and it rewrites 



Aa 



eSArea 



B 



E 



za^eB{n+ 1/2) 



cste 



(G3) 



n=0 2LI 



where we have introduced a cutoff Nc corresponding to 
the number of the Landau level for which continuum limit 
fails : za?eB{Nc + 1/2) ^ z, the width of the spectrum 
of the adjacency matrix. This gives Nc ^ \/{a'^eB) ^ 
4>o I {Be?). Using that the total number of bonds is related 
to the number of elementary plaquettes by S = \ x 

(# of plaquettes) we find Area = ^tan(7r/z) B. Finally 
the WL correction may be expressed with the Digamma 
functioiP : 



Act(6I< 1) 



a /TT 
— tan — 
27r \z 



(G4) 



In 



\Bd^ 



cste . 



where A/'[a;(T)] is the (algebraic) number of elementary 
plaquettes surrounded the trajectories. When one is in- 
terested in large time scale properties ^ a^), eq. (F21 1 
shows that the path integral over Brownian paths xir) 
in the square network can be replaced by a path integral 
over planar Brownian curves rij) as 



a;(T)^nctwork 

^a^tan(-) f I?f(r) e"' . (G6) 

2 \ Z / Jf(r)eplano 



9Af[x{T)] = 2eS^[a;(r)] coincides with the magnetic 
fiux enclosed by the trajectory multiplied by 2e, where 
^[a;(T)] is the algebraic area enclosed by the curve 
(.4[a;(r)] — a'^J\f[x{T)] for the square network). Finally 
we can rewrite ( G5 1 as 



Act(6'< 1) 



-a tan ( — 

. z 

r{t)—r 



V z 



(G7) 



/?(0)=F 

with D* — 1/2. Then we can use the well-known re- 
suliPSl ^ to get Aa{e < 1) ~ 5^ tan(f ) + 
SttBD'L'^ ) ~ (bo^)]'^^^^'^' which precisely coincides 
with (1^. 



3. Square network 

The WL correction can be expressed more explicitely 
for rational fiuxes. We recall briefiy a derivation due to 
Dougot & RammaPSl for the case of a square network of 
dimension V = x Ny. We start from 



Aaie) 



L., 



+ 4sinh(^/7a) 



coth(y7 a) 
1 



N^Ny 



Tr 



(G8) 



\Ar(7,0), 
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We label the vertices a = (n, m), n for the position along 
the horizontal direction and m along the vertical direc- 
tion. We can choose a Landau gauge : Oap = m6 for 
horizontal wires, a = (n,m) and /3 = (n + l,m), and 
Sap = for vertical wires. 

The computation of the MC for a square network 
submitted to a magnetic field requires to consider the 
problem of determination of the spectrum of Hap = 
—aai3e~^^°"^ (Hofstadter problem). In the Landau gauge 
where the flux is put along the n axis the solution of 
Hip = £tp can be written as V'n,™ — <&mG''°'"". We obtain 
the Harper equation <i>,„_i + [e + 2cos(fca; + ■md)]^m + 
= 0. 

If the flux is rational 6p,q — 2'np/q with p, g G N, the 
Harper equation is periodic with periodicity q. Therefore, 
writing m = r + qs with s € Z and r g {1, • • • , g} we can 
look for solutions of the form = (yS^e'*''''*. The wave 
function ipr is solution of a linear system Dq(p = where 
the hermitian q x g-matrix Dq is defined by 



m 



The secular equation det Dq = gives the q energy bands 
denoted £r{k). Interestingly the secular equation can 
be rewritten as det Dq = {—l)'^[Pp^q{e) — 2{cosqkx 
cosfcj,)] = 0, where Pp,q{£) is a polynomial of degrei 
q. For example = —e, -Pi, 2(e) = — 4, Pi, 3(e) 

—e^ + 6e, etc. 

The weak localization correction involves 



N,Ny 



— T f 1 
2. ^2^ , 1 



(27r)2 q ^ 4cosh(77a) + £^(fc) 



(GIO) 



2ir 



P;,(4cosh(V7a)) 



(27r)2 Pp_q(4cosh(y^a)) — 2cosfca. — 2cosfc, 



V 

(Gil) 



Dr+l,r 
I Dq^l 



e - 
1 



2 cos ( k- 



(fc. + 27r^) 



2nf] forr== I,--, 
for 7' = 1. 



,9-1 
(G9) 



1 



where we have used the equality Ylt=i 
for a polynomial of degree q : P(e) 
Integration over k leads to the result ( |75| of Dougot 
Rammal. 



" = ^^^^^ 

nLi(e.-^)- 
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